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This  study  Is  a  part  of  the  Advanced  Ballistics  Reentry 
Systems  (ABRES)  program,  a  major  concern  of  which  is  to  estimate 
nose  cone  erosion  by  wat?r  and  Ice  particles  In  natural  clouds. 

The  authors  acknowledge  with  gratitude  the  support  and 
guidance  of  Dr.  Robert  Cunningham  and  Mr.  Morton  Glass  of  the 
Convective  Cloud  Physics  Branch,  AFCRL. 
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I  NT  ROD  JCT I  ON 


It  is  current  practice  in  cloud  physics  research  to  make 
in  situ  measurements  of  hydrometeors  via  instrumented  low  speed  air¬ 
planes.  Quantitative  estimates  of  hydrometeor  concentrations  and  size 

spectra  are  obtained  with  particle  repl  icators^  and  Knollenberg 

(31 

particle  spectrometers'  .  Ideally,  these  Instruments  art  mounted  such 
that  measurements  are  made  (or  samples  taken)  in  the  undistrubed  free- 
stream,  thereby  avoiding  concentration  distortion  caused  by  airflow  about 
the  airplane.  However,  because  of  competition  for  the  best  locations  and 
the  many  problems  with  remote  control  of  instruments,  ideal  mounting  is 
frequently  not  possible.  Then  it  is  necessary  to  make  rational  choices 
between  the  various  options  available,  and  if  possible,  to  correct 
observed  data  to  remove  distortions.  A  method  that  has  proved  success¬ 
ful  in  accomplishing  both  of  these  goals  is  described  herein. 

The  method  has  been  applied  to  instrumentation  on  three  air¬ 
planes:  Lockheed  C130A  and  C130E  transports,  and  a  Cessna  Citation 
executive  jet.  The  Lockheed  airplanes  (Fig.  1),  which  are  similar  in 
appearance  except  for  a  longer  nose  radome  on  the  C130E,  are  instrumented 
for  use  by  the  Convect’ve  Cloud  Physics  Branch,  Cambridge  Research  Labora¬ 
tories.  We  have  done  extensive  theoretical  studies  of  formvar  particle 
replicators  mounted  on  these  airplanes.  The  Cessna  airplane  (Fig.  2) 
has  been  instrumented  for  cloud  physics  studies  by  Meteorology  Research, 
Inc.  We  have  studied  Knollenberg  particle  spectrometers  mounted  on  this 
airplane.  Flight  conditions  for  the  three  airplanes  are  given  in  Table  1. 


1.  P.  Spyers-Duran  and  R.  R.  Braham,  "An  Airborne  Continuous  Cloud 
Particle  Replicator,"  J.  Appl.  Meteor.  6,  1108  (1967). 

2.  J.  Hallett,  R.  W.  Hanaway,  and  P.  B.  Wagner,  "Design  and  Construc¬ 
tion  of  a  New  Cloud  Particle  Replicator  for  Use  on  a  Pressurized 
Aircraft,"  Desert  Research  Institute,  Reno,  Nevada,  AFCRL-72-0410 
(31  May  1972).  AD-753  091 . 

3.  R.  G.  Knollenberg,  "The  Optical  Array:  An  Alternative  to  Scattering 
or  Extinction  for  Airborne  Particle  Size  Determination,"  J.  Appl. 
Meteor.  9,  86  (1970). 
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FIGURE  2.  Cessna  Citation  Executive  Jet.  Wingspan  -  44  feet;  overall  length  -  44  feet; 
fuselage  radius  -  265  feet.  Location  of  the  Knollenberg  spectrometers  is 
shown  in  Fig.  25. 
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NATURE  OF  THE  PROBLEM 


Consider  axisymmetric  airflow  about  a  prolate  ellipsoid 
(Fig.  3),  which  we  can  use  as  a  simple  approximation  to  an  aircraft 
fuselage.  If  a  particle  sampling  instrument  were  positioned  near  the 
fuselage,  for  example  at  the  point  marked®  in  Fig.  3,  and  if  the 
particles  were  to  follow  the  stream  flow,  then  from  the  streamline 
spacing  contraction  it  is  apparent  that  a  higher  particle  flux  would  be 
measured  than  exists  in  the  free-stream.  This  could  be  evaluated  by 
straightforward  analysis  of  the  flow. 

The  situation  is  more  complicated  when  the  inertia  of  the 
particles  causes  them  to  deviate  from  the  streamlines.  Inertia  effects 
are  illustrated  by  Figs.  4,  5,  and  6,  which  show  computed  water  drop 
trajectories  for  drops  of  50,  100,  and  1000  ym  diameter,  respectively, 
in  the  potential  flow  field  of  Fig.  3.  In  Fig.  4  the  trajectories 
closest  to  the  ellipsoid  are  seen  to  deflect  outward.  In  Fig.  5  there 
is  both  impaction  on  the  ellipsoid  and  substantial  trajectory  deflection. 
These  effects  conspire  to  produce  quite  high  particle  fluxes  at,  for 
example,  point  1  in  Fig.  5,  whereas  at  point  2  we  should  not  expect  to 
sample  water  drops  of  this  size.  Point  2  lies  in  a  "shadow  zone"  region, 
and  adjacent  to  this  region  we  expect  high  particle  concentrations  and 
very  steep  concentration  gradients.  (Characteristically,  we  find  maximum 
concentration  distortion  for  water  drops  of  about  100  ym  diameter.)  For 
Fig.  6  we  see  that  considerably  larger  water  drops  possess  sufficient 
inertia  that  they  substantially  ignore  airflow  about  the  ellipsoid. 

Notice  the  left-to-right  downward  slant  of  the  trajectories  in 
Fig.  6.  This  is  caused  by  gravity  settling  of  the  drops.  Thus,  though 
the  airflow  is  axisymmetric,  the  particle  flux  is  not.  Moreover,  bodies 
of  revolution  are  poor  approximations  to  fuselages,  wings,  etc.,  and 
angle-of-attack  departures  from  ax  i  symmetric  flow  need  to  be  accounted 
for.  Clearly  three-dimensional  methods  are  required. 
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FIGURE  4.  Trajectories  of  50  um  diameter  water  drops  in  potential  airflow  about  a  prolate  ellipsoid  of 
fineness  ratio  2.  (The  ordinate  scale  is  expanded  by  a  factor  of  2.) 
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THEORY 


CONCENTRATION  FACTOR 

Principal  results  of  this  work  are  expressed  in  a  quantity 
called  concentration  factor.  Concentration  factor,  Cp  is  defined'4'  as 
the  ratio  of  particle  flux  (i.e.,  mass  of  particles  passing  per  second 
through  a  unit  area  normal  to  the  particle  velocity)  at  the  sampling  or 
target  point,  Ft,  to  the  particle  flux  in  the  free-stream,  F, 


0) 


The  ratio  of  particle  concentration  at  the  target  point  to  free-stream 
concentration,  C^,  Is 

CH  .  Cp  V/Vt  .  (2) 

where  V  is  free-stream  airspeed  ir-d  Vt  is  airspeed  at  the  target  point. 

In  this  latter  definition  we  ignore  difference  between  particle  and  air 
velocities.  Since  V^/V  =  1  (within  10X  for  all  cases  studied,  see 
Table  2)  we  confine  our  attention  hereafter  to  the  more  precise  quantity  Cp. 

In  three  dimensions  we  determine  concentration  factor  via  cal¬ 
culation  of  a  particle  flux  tube  (Fig.  7).  This  tube,  which  is  analogous 
to  a  streamtube,  is  determined  such  that  there  is  no  particle  flux 
through  its  boundaries;  therefore  mass  transfer  of  particles  is  equal 
through  all  cross-sections.  It  is  centered  about  a  central  trajectory  (the 
heavy  dashed  curve  in  Fig.  7)  that  passes  through  the  primary  target  point. 
The  initial  and  target  planes  are  perpendicular  to  the  central  trajectory. 


4.  R.  G.Dorsch  and  R.  J.  Brun,  "Variation  of  Local  Liquid-Water  Concen¬ 
tration  About  an  Ellipsoid  of  Fineness  Ratio  5  Moving  in  a  Droplet 
Field,"  NACA-TN-31 53  (July  1954). 
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TABLE  2 


RATIO  OF  AIRSPEED  AT  SAMPLING  POINTS  TO  FREE-STREAM  AIRSPEED 


Airplane 

Instrument 

Altitude 

(kft) 

Vt 

Lockheed  ) 

r  Particle 
replicator 

5 

1.102 

C130A  \ 

Particle 

replicator 

30 

1.102 

Lockheed 

C130E 

Particle 

replicator 

5 

1.040 

Cessna  , 

r  Precipitator 
particle 
spectrometer 

20 

1.017 

Citation 

}  Cloud 
particle 
^  spectrometer 

20 

1.020 

-  11  - 


If  M  is  the  particle  mass  transfer  rate  through  the  tube,  then 
at  any  point  along  the  tube 

M  =  FA  ,  (3) 

where  A  Is  the  perpendicular  cross  section  area  of  the  tube.  Since  M  Is 
constant  in  a  particle  flux  tube, 

CF  ”  ^ 

or 

Cp  =  lim  (t~\  »  (5) 

h  A->0\Mt/ 

At^0 

where  A  and  A{  are  the  cross-sectional  areas  of  the  flux  tube  in  the  free 
stream  and  at  the  target  point,  respectively. 

In  broad  outline,  our  procedure  is  as  follows.  We  define  a 
circular  "window"  (i.e.  flux  tube  surface  trace)  in  the  target  plane.  Our 
primary  target  point  is  at  the  center  of  this  circle.  A  number  of  evenly 
spaced  points  (usually  6  or  8)  are  chosen  on  the  window  circumference. 
Then,  by  use  of  an  iterative  procedure,  described  below  and  in  Appendix  A 
we  establish  the  particle  trajectories  that  pass  through  these  points. 

We  take  the  trajectory  intersection  points  with  the  target  and  initial 
planes  to  be  the  vertices  of  plane  polygons:  an  approximate  regular 
polygon  in  the  target  plane,  and  an  irregular  (distorted)  polygon  in  the 
initial  plane  (Fig.  8).  We  compute  the  areas  of  these  polygons,  and 
take  the  concentration  factor  to  be  the  ratio  of  these  areas. 

Owing  to  the  geometrical  connection  between  stream  tube  cross- 
sections  in  the  initial  and  target  planes,  we  hypothesize  that  they  can 
be  related  via  a  conformal  mapping  transformation.  Let  a  complex  point 


-  13  - 


Fuselage 
- #► 


Initial  plane 
wl ndow 


FIGURE  8.  Target  and  Initial  plane  flux  tube  cross  sections  for  50  ym 

water  drops  passing  through  the  C130A  particle  replicator  slit. 
Cp  =  1.271,  rw  -  0.02,  e  =  0.1 
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Zt  tar9et  plane  be  related  to  its  corresponding  point 

Z  (Z«c+U)  In  the  initial  plane  via  the  function  f(Zt).  That  Is 

Z  *  f(Zt)  .  (6) 

Our  objective  Is  to  find  the  po  nt  Z,  In  the  Initial  plane,  that  corre¬ 
sponds  to  a  given  target  point,  Zt  Q.  We  do  this  by  constructing  the  func¬ 
tion  f(Zt).  Since  f(Zt)  Is  an  inalytlc  function.  It  follows  that  a  Taylor's 
series  expansion  of  f(Zt)  about  the  target  point,  Zt  exists (  ,  viz. 


df(Z, 


f<zt>  ■  f<zt.o> +  <Vzt.0> 


<zt-zt.o>2,  t  ^<zt.o>  <V2t,/ 

,2  TT  ■  -»  ~  r-  *  ... 


dZ 


dZ*J 


m 


(7) 


We  begin  by  neglecting  all  terms  in  Eq.  (7)  higher  than  first 
order.  Then  we  have 


z  ■  8tC<Vzt.o> 


(8) 


where  B  and  C  are  complex  constants  and  Zfc  g  Is  the  given  point  in  the 
target  plane.  We  make  two  estimates  of  Z:  Z1  and  Z 2,  and  two  estimates 
of  the  corresponding  values  of  Zt:  Z  .  and  l.  This  yields  four 
simultaneous  equations: 


5.  G.  F.  Carrier,  M.  Krook,  C.  E.  Pearson,  Functions  of  a  Complex 
Variable  (McGraw-Hill  Book  Company,  196uTI  Sect.  2.6. 
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Z1  ■  6  ♦  c  <Zt,l  -  Zt.0> 


Z2 


B  ♦  C  <Zti2 


0) 


which  are  solved  for  B.  Our  next  estimate  of  Z  is  taken  to  be 


Z 


3 


=  B 


(10) 


By  use  of  ly  a  trajectory  is  calculated  and  Z^  3  is  found.  Equations  (9) 
are  solved  for  B(=Z^),  etc.  The  procedure  is  continued  until 


Here  e  is  a  prescribed  fractional  multiple  of  the  target  plane  window 
radius,  r^. 

This  start-up  procedure  uses  the  results  of  only  the  last  two 
trajectory  calculation'-  to  determine  the  next  set  of  Initial  coordinates. 
When  j=5,  a  least  squares  calculation,  which  uses  all  of  the  preceding 
results,  takes  over  and  is  used  thereafter.  Both  iteration  procedures 
are  described  in  more  detail  in  Appendix  A. 

PARTICLE  TRAJECTORY  CALCULATION 

The  equations  of  motion  of  a  heavy  particle  in  a  fluid  are 
based  on  the  assumption  that  the  bulk  fluid  flow  is  not  perturbed  by  the 
particles.  Therefore,  the  particles  move  under  influence  of  the  forces 
of  hydrodynamic  drag,  gravity,  buoyancy,  and  inertial  reaction  of  fluid 
carried  along.  For  particles  small  enough  for  application  of  Stokes 
drag  law,  the  theory  is  quite  adequately  developed.  For  larger  particles, 
the  theory  is  deficient  and  we  must  resort  to  approximate  methods^. 


6.  N.  A.  Fuchs,  The  Mechanics  of  Aerosols,  Translated  by  R.  E.  Daisley 
and  M.  Fuchs,  edited  by  C.  N.  Davies  (MacMillan,  New  York,  1964). 
Chapter  III. 
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Davies  and  Aylward^  computed  trajectories  of  small  particles 
In  flow  about  a  plate,  and  Davies  and  Peetz^  performed  similar  calcu¬ 
lations  for  flow  around  cylinders.  More  recently  accurate  calculations 
for  small  particle  flow  about  spheres^  and  plates^)  have  been 
reported . 

Calculations  for  spherical  particles  of  arbitrary  size  have 
been  reported  by:  Langmuir  and  Blodgett  for  flow  about  spheres,  cyl¬ 
inders,  and  ribbons^';  Dorsch  et  al .  for  flow  about  ellipsoids  of 

(121  fl3l 

revolution'  Norment  for  flow  In  nuclear  clouds'  Etkln  for  air 

1141  Ms) 

jets'  ' ;  Morsi  and  Alexander  for  flow  about  clylnders  and  airfoils' 

and  Chai  for  airfoils^).  All  of  these  calculations  are  two-dimensional. 


7.  C.  N.  Davies  and  M.  Aylward,  "The  Trajectories  of  Heavy,  Solid  Parti¬ 
cles  in  a  Two-Dimensional  Jet  of  Ideal  Fluid  Impinging  Normally  upon 
a  Plate,"  Proc.  Phys.  Soc.  (London)  B64,  889  (1951). 

8.  C.  N.  Davies  and  C.  V.  Peetz,  "Impingement  of  Particles  on  a  Trans¬ 
verse  Cylinder,"  Proc.  Roy.  Soc.  (London)  A234,  269  (1956). 

9.  K.  V.  Beard  and  S.  N.  Grover,  "Numerical  Collision  Efficiencies  for 
Small  Raindrops  Colliding  with  Micron  Size  Particles,"  J.  Atmos. 

Sci.  31_,  543  (1974). 

10.  R.  L.  Pitter  and  H.  R.  Pruppacher,  "A  Numerical  Investigation  of 
Collision  Efficiencies  of  Simple  Ice  Plates  Colliding  with  Super¬ 
cooled  Water  Drops,"  J  Atmos.  Sci.  31_,  551  (1974). 

11.  I.  Langmuir  and  K.  B.  Blodgett,  "Mathematical  Investigation  of  Water 
Droplet  Trajectories,"  General  Electric  Company,  Report  RL-225  (1945). 

12.  R.  G.  Dorsch,  R.  J.  Brun,  and  J.  L.  Gregg,  "Impingement  of  Water 
Droplets  on  an  Ellipsoid  with  Fineness  Ratio  5  in  Axisymmetric  Flow," 
NACA  TN-3099  (March  1954).  (Also  see  NACA  TN-2952,  2999,  3047, 

? ■ L5,  3410,  3153,  3586). 

13.  H.  G.  Norment,  "Research  on  Circulation  in  Nuclear  Clouds,  II,"  Technical 
Operations,  Inc., Report  TO-B  64-102  (November  1964).  AD-361  074. 

14.  B.  Etkin,  "Interaction  of  Precipitation  with  Complex  Flows,"  Proceedings 
of  the  Third  International  Conference  on  Wind  Effects  on  Buildings  and 
Structures  (Tokyo.  1971). 

15.  S.  A.  Morsi  and  A.  J.  Alexander,  "An  Investigation  of  Particle  Trajec¬ 
tories  in  Two-Phase  Flow  Systems,"  J.  Fluid  Mech.  55.,  193  (1972). 

16.  S.  Kuo-Kai  Chai,  "Droplet  Trajectories  Around  Aircraft  Wing,"  Thesis, 

U.  Nevada,  Reno  (November  1973). 
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If  the  particle  density  Is  large  compared  to  the  fluid,  which 
is  true  for  hydrometeors  In  air,  we  can  neglect  buoyancy  and  Inertial 
reaction  of  the  fluid  to  obtain  the  general  equation 

m  "  7  pAp  ■  V  |*f  -  »p|  cD  ♦  «*  02) 

where  m  Is  the  particle  mass,  A  the  particle  area  projected  In  the 

-►  P  -► 

direction  of  motion,  Vp  particle  velocity,  V^,  fluid  velocity,  CD  drag 
c 'efficient,  p  fluid  density,  and  g  gravity  acceleration.  Consider  a 
flow  of  constant  tree-stream  airspeed  V  around  a  body  of  characteristic 
dimension  L.*  T'ian  Eq.  (12)  can  be  non-dimensional  1  zed  to  yield 


■a? ~ "  (vfx  •  V> 

'  (vfy  '  V 

"  (Vfz  "  VDZ^  I  ;f  ■ 
a  tz  pz  TvTrN 


Here  length  is  scaled  by  L,  velocity  by  V,  and  time  by  L/V,  and 


(cdrn)/rn 

(14) 

V2/(Lg) 

(IS) 

*  Equivalent  results  are  obtained  by  assuming  either  a  moving  body  In  a 
stationary  fluid,  or  a  moving  fluid  about  a  stationary  body.  There¬ 
fore  we  use  whichever  concept  Is  most  expedient. 
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R M  -  Ei  v  -v.  V 
N  n  p  f 


Non-dimensional  quantities  are: 


V  vf 


particle  and  air  velocities 


CDRN  =  BN 


Froude  number 
Reynolds  number 
Best  number 

terminal  settling  speed  (Pj  is  computed  from  vT) 


Dimensioned  quantities  are: 

6  particle  dimension  (sphere  diameter  or  column 

base  width) 

p  air  density 

n  air  viscosity 

g  gravity  acceleration  constant 

V  free-stream  airspeed 

L  characteristic  dimension  of  body 

In  this  form,  the  equations  are  applicable  ta  any  flow  and  to  any  size 
and  shape  of  particle.* 


*  In  the  prior  work  cited  above,  riost  of  the  authors  use  some  modifica¬ 
tion  of  these  equations.  For  example,  even  though  their  trajectories 
are  essentially  horizontal  in  direction,  Dorsch  et  al.(]2).  Morsi  and 
AlexanderOB) f  and  Chaiv'B)  all  neglect  gravity.  Etkin'^J  uses  a 
constant  P.  Morsi  and  Alexander”*)  compute  and  use  P  independently 
for  each  space  component,  which  is  an  incorrect  procedure. 
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For  Stokes  drag  (R^<0. 1 )  P  has  the  constant  value  24.  For 
larger  R^,  P  Is  a  function  of  Reynolds  number  and  acceleration;  however 
the  dependence  on  acceleration  is  not  known.  It  Is  customary  practice 
to  use  steady-state  values  of  P,  which  are  determined  from  terminal 
settling  experiments.  Use  of  these  data  are  discussed  in  the  chapter  on 
Accuracy,  and  specifics  regarding  the  data  are  presented  In  Appendix  B. 

Equations  (13)  are  integrated  numerically  starting  at  a  point 
far  enough  upstream  that  essentially  free-stream  conditions  prevail. 

Krogh's  ordinary  differential  equation  Integrator  DVDQ^)  is  used.  This 
code  is  recomnended  by  Hull  et  al.^®^  and  gives  excellent  results  for 
this  problem.  The  technique  used  to  compute  vf  at  each  time  step  is 
described  below. 

THREE-DIMENSIONAL  FLOW  CALCULATION 

In  performing  concentration  factor  calculations  for  sampling 
sites  on  particular  airplanes  it  is  important  to  use  three-dimensional 
airflow.  This  is  the  only  way  to  adequately  account  for  particle  settling, 
airplane  geometry,  angle-of -attack,  airspeed  and  altlt  ide. 

Cloud  physics  airplanes  are  subsonic,  sampling  runs  being  made  , 
typically  between  100-150  kts.  indicated  airspeed.  Particle  measurement 
points  are  beyond  the  skin-friction  boundary  layer  (see  p.68),  and  should 
be  placed  to  avoid  separated  flow  regions.  Therefore,  potential  (i.e., 
frictionless,  incompressible,  laminar)  flow  calculations  are  quite 


17.  F.  T.  Krogh,  "Variable  Order  Integrators  for  Numerical  Solution  of 
Ordinary  Differential  Equations,"  Jet  Propulsion  Lab  Technology 
Utilization  Document  No.  CP-2308  (November  1970). 

18.  T.  E.  Hull,  W.  H.  Enright,  B.  M.  Fellen,  and  A.  E.  Sedgwick,  "Com¬ 
paring  Numerical  Methods  for  Ordinary  Differential  Equations," 

SIAM  J.  Numer.  Anal.  9,  603  (1972). 
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adequate.  We  use  a  code  developed  by  Hess  and  Smlth^®*^  for  calcula¬ 
ting  potential  flow  about  arbltr.  ry  three-dimensional  bodies.  (Recent, 
more  generalized  methods  reduce  .0  the  Hess-Smlth  procedure  for  compara¬ 
ble  application^^.)  The  Hess-imlth  code  requires  Input  of  a  digital 
description  of  the  aircraft  surface.  This  consists  of  the  coordinates 
of  the  corner  points  of  a  large  number  of  contiguous,  plane,  quadra- 
laterals.  An  example  of  the  digital  description  of  a  fuselage  section 
Is  shown  In  Fig.  9. 

The  Hess-Smlth  code  takes  each  quadralateral  panel  to  be  a 
uniform-distributed  source.  On  the  basis  of  the  boundary  condition 
that  there  be  zero  flux  through  the  centroid  of  each  panel,  and  given 
the  direction  of  the  free-stream  flow,  the  code  finds  the  source 
strengths  of  all  panels  by  Inversion  of  a  large  matrix  that  Includes 
all  possible  panel  interactions.  The  matrix  is  Inverted  only  once  for 
each  airplane  geometry,  provided  that  the  results  are  stored  for  future 
use. 

The  concentration  factor  calculations  require  flow  velocities, 
point-by-point  along  each  trajectory.  In  calculating  each  flow  velocity, 
contributions  from  all  panels  are  summed.  There  are  three  algorithms 
for  computing  contributions:  (1)  for  panels  that  are  close  to  the  calcu¬ 
lation  point,  a  detailed  calculation  is  used  that  accounts  for  exact 
panel  shape,  (2)  for  panels  at  intermediate  distances  a  multiple  expan¬ 
sion  is  used,  and  (3)  for  remote  panels  a  point  source  approximation  is 
used. 


19.  J.  L.  Hess  and  A.  M.  0.  Smith,  "Calculation  of  Non-Lifting  Poten¬ 
tial  Flow  about  Arbitrary  Three-Dimensional  Bodies,"  McDonnell 
Douglas  Report  E.  S.  40622  (15  March  1962).  AD-282  255. 

20.  J.  L.  Hess  and  A.  M.  0.  Smith,  'Calculation  of  Potential  Flow  About 
Arbitrary  Bodies,"  in  Progress  in  Aeronautical  Sciences,  Vol.  8, 
edited  by  D.  Kuchemann  (Perganroon  Press,  New  York,  1967). 

21.  F.  A.  Woodward,  "Analysis  and  Design  of  Wing-Body  Combinations  at 
Subsonic  and  Supersonic  Speeds,"  J.  Aircraft  5,  528  (1968). 
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FIGURE  9.  Computer-prepared  plot  of  the  digital  description  of  the 
nose  and  cabin  sections  of  the  Lockheed  C130A  airplane. 
®marks  the  location  of  the  particle  replicator. 
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To  perform  these  calculations,  we  have  developed  a  subroutine 
that  consists  of  various  extra  :ted  and  modified  portions  of  the  Hess- 
Smlth  code.  This  subroutine  1<.  generalized  such  that  given  the  source 
strength  results  for  any  three-dimensional  body,  It  will  provide  v^ 
for  any  Input  point  (x,y,z).  It  also  checks  each  Input  point  to  deter¬ 
mine  If  It  Is  Inside  the  airplane  body. 

Hess  and  Smith  ’•ave  essessed  the  accuracy  of  their  calcula¬ 
tions  for  a  number  of  body  shai  es  and  conditions;  excellent  agreement 
with  both  theory  and  experimen  is  found^’2^.  Our  assessment  is 
presented  In  the  chapter  on  Ac>  uracy.  Of  course  accuracy  depends 
on  the  fineness  of  resolution  (  f  the  panel  description  of  the  body. 

Here  some  compromise  is  called  for.  The  smaller  the  panels  the  finer 
the  resolution,  and  the  fewer  of  them  for  which  the  most  exacting  of 
the  three  algorithms  must  be  used.  On  the  other  hand,  the  number  of 
panels  increases  inversely  as  the  square  of  their  linear  size.  We  have 
used  che  following  criteria  in  setting  up  panel  structures:  For 
those  parts  of  the  airplane  traversed  by  particle  trajectories,  we  try 
to  keep  the  panel  edges  between  6"  to  8"  In  length  (sampling,  l.e.  target, 
points  are  9"  to  15"  from  the  fuselage).  Where  allowed  by  simplicity  of 
surface  shape,  remote  panels  can  be  larger.  Remote,  downstream  com¬ 
plexities,  such  as  the  wings  and  the  tail,  are  ignored.  The  cylindri¬ 
cal  portion  of  the  fuselage  is  extended  to  approximately  five  times 
the  length  of  the  nose  section,  as  recommended  by  Hess  and  Sm1thv“  ' . 

Computer  time  required  for  concentration  factor  calculation  is 
largely  dependent  on  the  number  of  v^  calculations  required.  On  the 
CDC  6600  computer,  one  calculation  requires  on  the  order  of  0.15  sec. 
The  number  of  vf  required  per  trajectory  varies  from  about  60  to  300.  A 
typical  number  of  trajectories  required  is  25.  Thus,  computing  time, 
even  on  a  large  computer,  can  be  considerable. 
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HYDROMETEOR  TYPES 


WATER  DROPS 

From  both  theoretical  and  experimental  viewpoints,  by  far  the 
simplest  hydrometeor  to  deal  with  Is  the  water  drop.  Accordingly,  con¬ 
fidence  In  concentration  factor  results  for  water  drops  Is  highest.  For 
this  reason  our  preliminary  studies  were  done  for  water  drops,  and  our 
"benchmark"  results  for  all  airplanes  are  for  water  drops. 

Terminal  settling  speeds  in  air  have  been  determined  experi- 
(22) 

mentally  by  Gunn  and  Kinzer'  and  confirmed  by  Beard  and 

(23) 

Pruppacher'  ' .  Water  drops  smaller  than  about  400  ym  In  diameter  are 

(24  251 

essentially  spherical'  *  ' .  Since  concentration  factors  for  water 

drops  larger  than  400  ym  are  close  to  unity,  we  use  drag  data  for  spheres 
in  our  trajectory  calculations.  In  this  way  we  avoid  having  to  cope 
with  drop  deformation  effects  at  high  altitudes,  for  which  little  data 
exists. 

ICE  CRYSTALS 

From  the  ABRES  Program  viewpoint,  interest  in  ice  hydrometeors 
is  probably  greatest.  Within  limits  imposed  by  local  cloud  conditions, 


22.  R.  Gunn  and  G.  D.  Kinzer,  "The  Terminal  Velocity  of  Fall  for  Water 
Droplets  in  Stagnant  Air,"  J.  Meteor.  6.  243  (1949). 

23.  K.  V.  Beard  and  H.  R.  Pruppacher,  "A  Determination  of  the  Terminal 
Velocity  and  Drag  of  Small  Water  Drops  by  Means  of  a  Wind  Tunnel," 

J.  Atmos.  Sci.  26,  1066  (196  9). 

24.  H.  R.  Pruppacher  and  K.  V.  Beard,  "A  Wind  Tunnel  Investigation  of 
the  Internal  Circulation  and  Shape  of  Water  Drops  Falling  at  Terminal 
Velocity  in  Air,"  Quart.  J.  Roy.  Met.  Soc.  96,  247  (1970). 

25.  H.  R.  Pruppacher  and  R.  L.  Pitter,  "A  Semi -Empirical  Determination 
of  the  Shape  of  Cloud  and  Rain  Drops,"  J.  Atmos.  Sci.  28,  86  (1971). 
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(26) 

Ice  hydrometeors  can  assume  any  of  a  vast  variety  of  geometrical  forms v  '. 
Nevertheless,  certain  systematlcs  and  trends  have  been  established,  and 
these  can  be  used  to  reduce  this  problem  to  manageable  portions.  For 
example,  under  natural  conditions  water  freezes  to  form  crystals  with 
hexagonal  symmetry.  Crystallization  habits  tend  toward  hexagonal  based 
columns  or  plates.  Moreover,  whether  plates  or  columns  are  produced 
seems  to  depend  mostly  on  temperature^26"29^.  More  complex  ice  hydro¬ 
meteor  forms  are  either  variations  of  these  simple  forms,  more  complex 
combinations  of  them,  or  very  complex  and  essentially  Infinitely  variable, 
dendritic  forms. 

To  date,  our  work  has  concentrated  on  columns.  There  is  con¬ 
siderable  interest  In  this  crystal  form,  especially  since  they  are  found 

•  ( 30) 

to  be  the  major  constituent  of  cirrus  clouds'  . 

A  number  of  studies  of  the  dimensions  of  columnar  ice  have  been 

(31 ) 

reported.  The  most  complete  is  that  of  Auer  and  Veal'  ' ,  and  we  have 
chosen  to  use  their  results.  Their  equations  relating  the  length,  i,  and 
base  width,  6,  (see  Fig.  10)  of  natural  ice  columns  are  (dimensions  in  urn): 


26.  C.  Magono  and  C.  W.  Lee,  J.  Fac.  of  Sci., Hokkaido  U.,  Ser.  VII, 

Vol .  II,  321  (1966). 

27.  N.  H.  Fletcher,  The  Physics  of  Rainclouds,  (Cambridge  University 
Press,  1966). 

28.  B.  J.  Mason,  The  Physics  of  Clouds  (Clarendon  Press,  Oxford,  1971). 

29.  A.  Ono,  "Growth  Mode  of  Ice  Crystals  in  Natural  Clouds,"  J.  Atmos. 
Sci. 27,  649  (1970). 

30.  A.  J.  Heymsfield  and  R.  G.  Knollenberg,  "Properties  of  Cirrus 
Generating  Cells,"  J.  Atmos.  Sci.  29,  1358  (1972). 

31.  A.  H.  Auer  and  D.  L.  Veal,  'The  Dimension  of  Ice  Crystals  in 
Natural  Clouds,"  J.  Atmos.  Sci.  27,  919  (1970). 
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End-on 


Ax is -on 

drag 

force 


Volume: 

Vol  =  6zi  =  0.64951 96 2£. 


Mass: 

m  =  0.64951 962£  pp 

Diameter  of  Water  Drop  of  Equal  Mass: 

4„  =  1 .0744786(4*4  Pp/Pwater) 

FIGURE  10.  Properties  of  hexagonal -based  plates  and  columns 
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6  =  -  8.479  +  1.002*  -  0.00234a2  ;  A  <  200  pm 

6  =  11. 3A0*414  ;  i  >  200  pm  (17) 

Densities  for  solid  and  hollow  ice  columns  were  measured  by 
(321 

Jayaweera  and  Ryan'  ' .  They  computed  volumes  from  the  dimensions  of 
the  individual  crystals,  and  masses  were  estimated  from  the  melted 

3 

water  drop  diameters.  The  resulting  densities  were  0.7  g/cm  for  solid 

3 

columns,  and  0.36  g/cm  for  hollow  columns.  These  densities  are  used 
in  our  calculations. 

In  free-fall  settling  experiments,  columns  are  observed  to 
orient  their  axes  horizontally^22^ .  Moreover,  Bragg  et  al.^24^  find 
for  Reynolds  numbers  greater  than  about  .05,  that  this  orientation  is 
assumed  very  rapidly.  Accordingly,  we  take  our  column  settling  speed 
to  be  that  for  the  horizontal  orientation. 

During  trajectory  calculations,  we  take  the  drag  force  vector 
to  be  parallel  with  vf-vp.  Drag  coefficient  data  exist  for  two  column 
orientations  relative  to  the  drag  force  vector  (see  Appendix  B).  These 
orientations  are  shown  in  Fig.  10.  In  light  of  the  overwhelming  ten¬ 
dency  of  columns  to  fall  with  their  axes  horizontal,  it  is  apparent 
that  the  "axis-on"  orientation  is  the  significant  one.  Nevertheless, 
we  have  computed  concentration  factors  for  both  orientations. 

It  is  important  to  realize  that  exceptions  to  the  above  are 
common.  Substantial  departures  from  stable  horizontal  settling  are 


32.  K.  O.  L.  F.  Jayaweera  and  B.  F.  Ryan,  "Terminal  Velocities  of  Ice 
Crystals,"  Quart.  J.  Roy.  Met.  Soc.  98,  193  (1972). 

33.  K.  0.  L.  F.  Jayaweera  and  B.  J.  Mason,  "The  Behavior  of  Freely 
Falling  Cylinders  and  Cones  in  a  Viscous  Fluid,"  J.  Fluid  Mech.  22, 
709  (1965). 

34.  G.  M.  Bragg,  L.  van  Zuider,  and  C.  E.  Hermance,  "The  Free-Fall  of 
Cylinders  at  Intermediate  Reynolds  Numbers,"  Atmos.  Environ.  8, 

755  (1974). 
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reported  by  Jayaweera  and  Mason  for  asymmetrically 
and  by  Zlkmunda  and  Vall^^^  for  rimed  Ice  columns 
field. 


loaded  model s^^, 
observed  In  the 


At  any  particular  point  In  space,  such  as  a  measurement 
(target)  point.  It  is  not  possible  to  predict  with  assurance  the 
orientation  of  a  column  axis  In  the  horizontal  plane.  Our  best  assump¬ 
tion  Is  that  this  orientation  Is  variable.  A  Knollenberg  particle 
spectrometer  will  register  a  projection  of  the  column  dimension  in  the 
direction  perpendicular  to  Its  linear  sensing  array.  If  we  assume  a 
random  column  orientation  In  the  horizontal,  then  we  show  In  Appendix  C 
that  the  ensemble  mean  projected  dimension  of  a  column,  <v>,  is 


<v>  =  j  (i  +  6) 


(18) 


where  i  and  6  are  as  defined  In  Fig.  10.  We  feel  that  <v>  provides 
reasonable  approximations  to  mean  values  of  column  dimensions 
registered  by  a  Knollenberg  spectrometer. 

Properties  of  the  Ice  columns  studied  are  given  in  Table  3. 
Work  on  Ice  plates  Is  underway,  but  not  far  enough  along  to  allow  us 
to  Include  results  here.  We  also  plan  to  treat  other,  more  complicated 
shapes  such  as  spatial  dendrites. 


35.  K.  0.  L.  F.  Jayaweera  and  B.  J.  Mason,  "The  Falling  Motions  of 
Loaded  Cylinders  and  Discs  Simulating  Snow  Crystals,"  Quart.  J. 
Roy.  Met.  Soc.  92,  151  (1966). 

36.  J.  Zikmunda  and  G.  Vali,  "Fall  Patterns  and  Fall  Velocities 
of  Rimed  Ice  Crystals,"  J.  Atmos.  Scl.  29,  1334  (1972). 
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COLUMNAR  CRYSTAL  PROPERTIES 
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Solid  column  density  is  taken  to  be  0.7  g/cm  ,  and  hollow  column  density  is  taken  to  be  0.36  g/cm  . 


RESULTS 


PARTICLE  REPLICATOR  ON  THE  LOCKHEED  C130A 

(21 

The  particle  replicator  arm'  '  exits  the  C130A  fuselage  in  the 
cabin  section  just  aft  of  the  cockpit,  at  the  location  marked  in  Fig.  9 
and  as  shown  in  Fig.  11.  It  is  mounted  perpendicular  to  the  fuselage 
symmetry  plane.  As  determined  by  measurement  on  the  airplane,  the  parti¬ 
cle  intake  slit  is  14.3  inches  from  the  fuselage,  measured  along  the  arm. 

Our  digital  description  of  the  C130A  fuselage  was  developed 
from  Lockheed  engineering  drawings  of  moldline  contours.  Accuracy  of 
individual  coordinates  is  several  hundredths  of  an  inch.  Aft  of  the 
cabin  section,  the  fuselage  is  formed  from  sections  of  two  circular 
cylinders.  The  complete  fuselage  description  is  shown  in  Fig.  12.  A 
tot.1  of  1264  quadralateral  panels  are  used  in  this  description. 

Engineering  drawings  of  airplanes  are  based  on  a  coordinate 
system  that  appears  to  be,  in  part,  a  relic  of  shipbuilding.  The  fuse¬ 
lage  axis  coordinate  is  called  the  fuselage  station,  FS;  perpendicular 
to  the  FS  axis  we  have:  in  the  horizontal  the  buntline  axis,  BL,  and 
in  the  vertical  the  waterlevel  axis,  WL.  The  C130A  nose  tip  is  at 
FS=61",  the  bulkhead  that  separates  the  cabin  from  the  cylindrical  part 
of  the  fuselage  is  at  FS=245",  the  fuselage  synmetry  plane  is  at  BL=0, 
and  the  cabin  flight  deck,  as  well  as  the  center  of  the  principal  circu¬ 
lar  cylinder  (of  radius  85")  is  at  WL=200".  In  this  system  the  replica¬ 
tor  slit  has  coordinates: 

FS  =  17f .87" 

BL  =  85.45" 

WL  =  225.74" 

In  the  airflow  and  trajectory  computations,  we  take  the  origin  at: 

FS  =  245" 

BL  =  0 
WL  =  200" 
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FIGURE  11.  Lockheed  C130A  with  formvar  replicator  arm 
in  position  (flowplate  not  mounted) 
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The  characteristic  dimension  of  the  airplane,  L  (see  discussion  of 
Eqs.  (13)),  is  taken  as  the  length  of  the  cabin  section,  L=184".  This 
is  the  distance  from  the  tip  of  the  nose  to  the  cabin  bulkhead.  All 
linear  measures  are  normalized  by  this  length. 

Concentration  factors  for  water  drops  at  and  near  the  replica¬ 
tor  slit  at  5  kft,  and  at  the  slit  at  30  kft,  are  shown  In  Fig.  13  and 
listed  in  Table  4.  Flight  conditions  are  given  In  Table  1.  According 
to  these  results  the  replicator  is  mounted  In  a  relatively  favorable 
place,  In  the  sense  that  drops  of  all  sizes  are  sampled;  but  certainly 
it  is  not  free  of  concentration  distortion.  This  distortion  is  clearly 
evident  in  the  stereographic  plots  shown  in  Fig.  14.  In  these  plots  an 
enlarged  particle  flux  tube  is  shown  as  It  reaches  the  replicator  slit. 
Contraction  of  the  tube  relative  to  its  upstream  shape  is  apparent  (also 
see  Fig.  8).  Measured  fluxes  are  about  4G%  higher  than  their  free-stream 
values  for  drops  between  60-100  ym  diameter,  depending  on  altitude.  This 
represents  concentration  distortion  of  about  30%  as  obtained  by  dividing 
Cp  by  vT  (Table  2).  Moreover,  )nly  very  small  and  very  large  drops  are 
free  of  the  effect.  The  shift  of  concentration  factor  peak  to  smaller 
particle  sizes  as  altitude  is  increased  is  a  typical  trend. 

Concentration  factors  for  solid  and  hollow  ice  columns  at  the 
replicator  slit  at  5  kft  altitude  are  plotted  in  Figs.  15-18  and  listed 
in  Tables  5  and  6.  Properties  of  the  ice  crystals  are  given  in  Table  3. 
Calculations  were  done  for  both  the  "axis-on"  and  "end-on"  orientations 
with  respect  to  the  drag  vector  (Fig.  10).  Terminal  settling  with  long 
axis  horizontal  was  used  for  both  orientations.  As  is  evident  in  Figs. 
15-18,  there  are  large  differences  in  results  for  the  two  orientations, 
the  end-on  concentration  factors  for  the  large  columns  being  considerably 
smaller.  As  noted  above  (p.  27),  the  axis-on  orientation  is  observed  to 
be  the  stable  one  in  free-fall  settling  experiments.  Moreover,  Bragg 
et  al.^)  observe  that  this  orientation  is  assumed  very  rapidly  for 
columns  large  enough  to  be  of  interest  here.  Therefore,  we  consider  only 
the  axis-on  results  to  be  significant,  and  ignore  the  end-on  results 
hereafter. 
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CONCENTRATION  FACTOR 


FIGURE  13.  Water  drop  concentration  factors  at  and  near  the  particle 
replicator  slit  on  the  Lockheed  C130A 
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TABLE  4 


WATER  DROP  CONCENTRATION  FACTORS  AT  AND  NEAR  THE 
REPLICATOR  SLIT  ON  THE  LOCKHEED  C130A 
(See  Table  1  for  flight  conditions) 


Concentration  Factors _ 

30  kft 

5  kft  altitude _  altitude 


Drop 

diameter 

M 

Slit 

10  cm 
inboard* 
from  slit 

10  cm 
outboard* 
from  slit 

Slit 

10 

1.132 

15 

1.136 

20 

1.153 

1.184 

1.126 

1.176 

25 

1.206 

30 

1.186 

.241 

1.153 

1.236 

40 

1.226 

1.317 

1.180 

1.300 

50 

1.271 

1.412 

1.200 

1.346 

60 

1.225 

1.365 

70 

1.346 

1.600 

1.253 

1.353 

85 

1.389 

1.677 

1.255 

1.314 

100 

1.397 

1.693 

1.266 

1.277 

120 

1.623 

125 

1.370 

1.223 

150 

1.327 

1.509 

1.230 

1.179 

200 

1.252 

1.358 

1.192 

1.121 

300 

1.162 

1.208 

1.135 

1.070 

500 

1.072 

1.093 

1.058 

1.024 

1000 

1.010 

.998 

*  Distances  from  the  replicator  slit  are  along  the  replicator 
arm. 
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FIGURE  14a.  Stereographic  plots  of  an  eight-trajectory,  20  ym-diameter  water  drop  flux  tube 
to  the  Lockheed  C130A  particle  replicator  slit.  The  central  trajectory  is  also 
shown.  Cr-  =  1.15. 


Stereographic  plots  of  a  six-trajectory,  100  pm-diameter  water  drop  flux  tube  to 
the  Lockheed  C130A  particle  replicator  slit.  The  central  trajectory  also  is  shown. 

r  =  .02,.  e  =  .05,  Cc  =  1.41. 
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FIGURE  15.  Concentration  factor  v 
replicator  slit  on  the 
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MEAN  PROJECTED  COLUMN  DIMENSION,  <v>,  (ym) 

FIGURE  17.  Concentration  factor  vs.  mean  projected  column  dimension  for  solid  ice  columns  (p  =  0.7  g/cm3) 
at  the  particle  replicator  slit  on  the  Lockheed  C130A.  Altitude  is  5  kft.  p 
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FIGURE  18.  Cone 


TABLE  5 


ICE  COLUMN  CONCENTRATION  FACTORS  AT  THE  PARTICLE  REPLICATOR 
SLIT  ON  THE  LOCKHEED  C130A  (AXIS-ON  DRAG) 

Altitude  ■  5  kft 
(see  Tables  1  and  3) 


Length,  l 
(vm) 

Width,  6 
(uni) 

Concentration 

Factors 

Solid  Columns 

Hollow  Columns 

50 

35.8 

1.207 

70 

50.2 

1.260 

1.199 

85 

59.8 

1.291 

1.223 

100 

68.3 

1.323 

1.244 

200 

98.3 

1.416 

1.329 

300 

117.0 

1.451 

1.375 

500 

144.0 

1.443 

1.426 

700 

166.0 

1.423 

1.439 

900 

184.0 

1.408 

1.467 

1500 

227.0 

1.348 

1.453 

TABLE  6 


ICE  COLUMN  CONCENTRATION  FACTORS  AT  THE  PARTICLE  REPLICATOR 
SLIT  ON  THE  LOCKHEED  CI30A  (END-ON  DRAG) 

Altitude  ■  5  kft 
(see  Tables  1  and  3) 


Length,  fc 
(vm) 

Width,  6 
(ym) 

Concentration 

Factors 

Solid  Column 

Hollow  Column 

50 

35.8 

1.183 

1.151 

70 

50.2 

1.244 

1.185 

85 

59.8 

1.276 

1.212 

100 

68.3 

1.295 

1.235 

120 

78.1 

1.278 

1.264 

140 

85.9 

1.243 

1.282 

160 

91.9 

1 .206 

1.284 

200 

98.4 

1.155 

1.272 

300 

117.0 

1.151 

1.274 

380 

129.1 

1 .142 

1.269 

-  43  - 


Water  drop  results  « 1  so  are  plotted  In  Figs.  15-18  for  compari¬ 
son  with  those  of  the  Ice  columns.  We  see  that  the  ice  curves  peak  at 
larger  Cp  values  and  for  larger  particles  than  for  water  drops.  In  terms 
of  the  mean  projected  column  dimension  (see  p.  28).  the  curves  peak  at 
^300  pm  for  solid  columns  and  a-800  ym  for  hollow  columns.  Maximum  con¬ 
centration  distortion  of  about  30%  Is  Indicated,  with  slightly  smaller 
distortions  distributed  over  a  wide  range  of  particle  sizes. 

PARTICLE  REPLICATOR  ON  THE  LOCKHEED  C130E 

On  this  airplane  the  replicator  arm  exits  the  fuselage  through 
a  window  in  the  forward  cargo  compartment.  The  replicator  arm  is  perpen¬ 
dicular  to  the  fuselage  axis,  but  not  horizontal;  its  geometry  is  shown 
in  Fig.  19.  The  coordinate  system  (FS,  BL,  WL)  is  identical  to  the  one 
used  for  the  C130A,  and  we  use  the  same  origin  and  scaling  factor  for  our 
computations  (see  pp.  30,  33).  The  replicator  location  is  shown  in  Fig.  20. 

The  C130A  digital  description  was  modified  to  allow  for  the 
C130E  nose  radome  and  the  different  trajectory  paths.  Quadralateral 
spacing  is  as  discussed  for  the  C130A.  A  total  or  1692  quadralateral s 
are  used  in  this  description.  The  complete  fuselage  is  shown  in  Fig.  21. 

Concentration  factor  calculations  were  made  for  water  drops  at 
5  kft  altitude.  Flight  conditions  are  given  in  Table  1.  Results  are 
shown  in  Fig.  22  and  listed  in  Table  7.  Our  calculations  indicate  that 
water  drops  over  a  diameter  range  of  about  90-300  ym  cannot  be  sampled 
at  the  replicator  slit.  The  calculations  crudely  define  the  shadow  zone 
indicated  by  the  shaded  region  in  Fig.  22.  The  replicator  arm  is  being 
lengthened  by  8"  to  alleviate  this  problem;  however,  even  at  this  distance 
the  calculations  indicate  substantial  concentration  distortion. 

Stereographic  trajectory  plots  are  shown  in  Fig.  23.  These 
trajectories  are  to  a  point  10  cm  outboard  from  replicator  slit,  measured 
along  the  replicator  arm.  Deflection  and  flux  distortion  are  evident, 
the  concentration  factor  at  this  point  being  2.5.  The  cause  of  the 
shadow  zone  is  obvious.  Trajectories  closest  to  the  airplane  almost 
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DISTANCE  FROM  THE  FUSELAGE  (inches) 


WATER  DROP  DIAMETER  (vm) 


FIGURE  22.  Concentration  factor  contours  vs.  water  drop  diameter 
along  the  particle  replicator  arm  on  the  Lockheed  C130E 


FI6URE  23a.  Stereographic  plots  of  a  six-trajectory,  150  um-diameter  flux  tube  to  a  point  10  cm 
outboard  from  the  Lockheed  C130E  particle  replicator  slit.  The  central  trajectory 
also  is  shown. 
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FIGURE  23b.  Stereographic  plots  of  a  six-trajectory,  150  Mm-diameter  flux  tube  to  a  point  10  cm 
outboard  from  the  Lockheed  C130E  particle  replicator  slit.  The  central  trajectory 
also  is  shown. 


TABLE  7 


WATER  DROP  CONCENTRATION  FACTORS  FOR  THE  REPLICATOR 
ON  THE  LOCKHEED  C130E 

(see  Table  1  for  flight  conditions) 


Water 

drop 

diameter 

Distance  from  the  Fuselage  Along  the 
Replicator  Arm  (Inches) 

Slit 

haul _ 

13.75 

17.687 

21.624 

25.561 

29.489 

50 

1.265 

1.187 

1.148 

1.123 

1.106 

70 

^1 .40 

1.319 

1.228 

1.176 

1.145 

100 

— 

1.630 

1.366 

1.262 

1.204 

150 

— 

2.460 

1.568 

1.370 

1.275 

200 

— 

2.299 

1.578 

1.386 

1.289 

300 

-- 

1.719 

1.443 

1.322 

1.248 

500 

1.401 

1.272 

1.202 

1.154 

1.117 

600 

1.272 

700 

1.201 

800 

1.157 

1.1  7 

1.089 

1.069 

1.052 

-  51  - 


graze  the  fuselage.  It  Is  clear  that  drops  of  this  size  cannot  reach 
the  replicator  since  they  are  intercepted  by  the  fuselage  in  an  area 
below  the  cabin. 

KNOLLENBERG  PARTICLE  SPECTROMETERS  ON  THE  CESSNA  CITATION 

The  Knollenberg  particle  spectrometers  are  mounted  on  the 
emergency  exit  door  In  the  forward  part  of  the  passenger  cabin  (Figs.  24 
and  25).  The  precipitation  particle  spectrometer  has  a  particle  size 
range  200-3000  um,  and  the  cloud  particle  spectrometer  has  a  particle 
size  range  20-300  ym.  Measurement  points  are  9"  from  the  fuselage  for 
both  instruments.  Both  instruments  are  mounted  with  their  axes  perpen¬ 
dicular  to  the  fuselage  symmetry  plane. 

The  fuselage  description  (Figs.  25  and  26)  was  developed  from 
dimensioned,  line  drawings  taken  from  a  brochure.  The  drawings  were 
enlarged  and  copied  to  graph  paper  by  use  of  a  Kargel  Reflecting  Pro¬ 
jector.  Individual  coordinate  values  are  accurate  to  within  2V, 
though  rather  extensive  interpolation  was  required  in  some  sections. 

On  most  of  the  fuselage,  especially  the  parts  traversed  by  particles, 
quadralateral  panel  edges  were  kept  to  about  6".  A  total  of  1304  panels 
are  used  in  the  complete  fuselage  description. 

The  nose  and  cockpit  part  of  the  fuselage  is  10'  long.  The 
passenger  cabin  is  a  circular  tube  with  a  radius  of  2.65'.  The  particle 
spectrometers  are  14'  aft  of  the  nose  tip. 

Coordinate  measurements  were  taken  in  a  (FS,  BL,  WL)  coordin¬ 
ate  system  with  FS=0  at  the  nose  tip,  and  the  fuselage  axis  at  BL=0, 
WL=4.25'.  In  this  system  the  centers  of  the  optical  paths  of  the  parti¬ 
cle  spectrometers  are: 

precipitation  particle  spectrometer  - 
FS  =  13.95' 

BL  =  3.40' 

WL  =  4.25' 
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FIGURE  24.  Cloud  physics  instrumentation  mounted  on  the  Cessna  Citation. 

The  uppermost  and  lowermost  instruments  are  the  cloud  particle 
and  precipitation  particle  spectrometers,  respectively. 

(Photo  courtesy  of  Meteorology  Research,  Inc.) 


FIGURE  25.  Computer-prepared  plot  of  the  digital  description  of 
the  forward  fuselage  of  the  Cessna  Citation.  The 
precipitation  particle  spectrometer  location  is  marked 
with#  and  the  cloud  particle  spectrometer  with  ff 
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cloud  particle  spectrometer  - 
FS  =  13.95' 

BL  =  3.056' 

WL  =  5.807' 

For  computer  calculations,  the  origin  of  coordinates  Is: 

FS  =  10' 

BL  =  0 
WL  =  4.25' 

The  characteristic  length  of  the  airplane,  L,  by  which  all  linear 
measures  are  normalized  (see  discussion  of  Eqs.  (13)),  is  10  feet. 

Water  drop  concentration  factors  for  the  particle  spectro¬ 
meters  are  plotted  in  Fig.  27  and  listed  in  Table  8.  Flight  conditions 
are  given  in  Table  1.  Extensive  calculations  also  were  done  for  ice 
columns.  Results  are  plotted  in  Figs.  28-31  and  are  listed  in  Table  9. 

Serious  problems  exist  at  both  spectrometers,  but  particularly 
at  the  cloud  particle  spectrometers.  From  Fig.  32  we  see  that  particles 
must  undergo  considerable  trajectory  deflection  to  reach  the  target 
points.  At  the  precipitation  particle  spectrometer,  concentration  dis¬ 
tortion  ranges  from  60-80*.  For  the  cloud  particle  spectrometer,  >ue 
deflection  point  is  closer  to  the  spectrometer  and  the  deflection  is  more 
abrupt,  which  explains  the  more  serious  problem  at  this  spectrometer.  At 
the  cloud  particle  spectrometer  a  narrow  shadow  zone  between  00- 1 20  pm 

is  indicated  for  water  drops.  For  ice  columns  at  the  cloud  particle 
spectrometer,  there  is  a  shadow  zone  above  <300  pm  mean  projected  dimen¬ 
sion  for  solid  columns,  and  ^00  pm  for  hollow  columns.  Large  concen¬ 
tration  distortions  are  indicated  for  smaller  columns. 

The  situation  here  is  similar  to  that  of  the  C130E  replicator. 
The  target  points  are  far  removed  from  the  free  stream:  being  too  far 
aft  and  too  close  to  the  fuselage.  Figure  32  suggests  the  possibility 
of  a  flow  separation  problem  at  the  cloud  particle  spectrometer, 
though  our  method  does  not  allow  Investigation  of  this  possibility.  It 
is  clear  that  both  spectrometers,  particularly  the  cloud  particle  spec¬ 
trometer,  are  very  unfavorably  mounted. 
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CONCENTRATION  FACTOR 


WATER  DROP  DIAMETER  (pm) 

FIGURE  27.  Concentration  factors  vs.  water  drop  diameters  for  the  Knollenberg 
particle  spectrometers  on  the  Cessna  Citation.  The  point  marked',* 
Is  hand-calculated  from  a  partial  result. 
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TABLE  8 


WATER  DROP  CONCENTRATION  FACTORS  FOR  THE  PARTICLE 
SPECTROMETERS  ON  THE  CESSNA  CITATION 

(see  Table  1  for  flight  conditions) 


Water 

drop 

diameter 

Precipitation 

Particle 

Spectrometer 

Cloud 

Particle 

Spectrometer 

30 

1.079 

1.126 

50 

1.198 

1.328 

60 

1.281 

70 

1.380 

1.712 

80 

1.480 

2.076 

90 

1.564 

95 

2.827 

100 

1.622 

0,5.6* 

no 

1.645 

120 

1.614 

150 

1.482 

2.286 

170 

1.880 

200 

1.324 

1.601 

300 

1.182 

1.270 

400 

1.112 

1.159 

600 

1.046 

1.066 

800 

1.018 

*  Estimated  from  partial  results 
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FIGURE  29.  Concentration  factor  vs.  particle  mass  for  solid  and  hollow  ice  columns  at  the 
cloud  particle  spectrometer  on  the  Cessna  Citation 


Hollow  ice  columns 


MEAN  PROJECTED  COLUMN  DIMENSION,  <v>,  (um) 

FIGURE  30.  Concentration  factor  vs.  mean  projected  column  dimension  for  solid  and  hollow  ice 
columns  at  the  precipitation  particle  spectrometer  on  the  Cessna  Citation 


CONCENTRATION  FACTOR 


FIGURE  31.  Concentration  factor  vs.  mean  projected  column  dimension 
for  solid  and  hollow  ice  columns  at  the  cloud  particle 
spectrometer  on  the  Cessna  Citation. 
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TABLE  9 


CONCENTRATION  FACTORS  FOR  ICE  COLUMNS  AT  THE  KNOLLENBERG 
SPECTROMETERS  ON  THE  CESSNA  CITATION 
(see  Tables  1  and  3) 


Length,  l 
(pm) 

Width,  6 
(u*n) 

Concentration 

Factors 

Precipitation 

Particle 

Spectrometer 

Hollow  Solid 

Cloud  Particle 
Spectrometer 

Hollow  Soli 

50 

35.75 

1.09 

1.08 

1.13 

70 

50.19 

1.15 

1 .12 

1.23 

85 

59.76 

1.21 

1.15 

100 

68.3 

1.25 

1.19 

1.41 

200 

98.4 

1.24 

1.51 

1.34 

1.90 

300 

117.0 

1.32 

1.66 

1.48 

2.48 

500 

144.5 

1.44 

1.76 

1.78 

700 

165.9 

1.76 

2.09 

900 

183.6 

1.65 

1.70 

2.36 

1000 

197.3 

1.75 

1 .68 

2.65 

1500 

226.5 

1.84 

1 .55 

2000 

262.8 

1.81 

1.43 

3000 

310.9 

1.71 

1.33 

4000 

223.1 

1.31 
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FIGURE  32a.  Stereographic  plots  of  four-trajectory,  150  yin-diameter  water  drop  flux  tubes  to 

the  Cessna  Citation  cloud  particle  spectrometer  (upper,  Cp  =  2.29)  and  precipitation 
particle  spectrometer  (lower,  Cp  =  1.48).  Individual  trajectories  are  not  resolved. 
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FIGURE  32b.  Stereographic  plots  of  four-trajectory,  150  ym-diameter  water  droD  flux  tubes  to 

the  Cessna  Citation  cloud  particle  spectrometer  (upper,  Cp  =  2.29;  and  precipitation 
particle  spectrometer  (lower,  Cp  =  1.48).  Individual  trajectories  are  not  resolved. 


ACCURAC  t  OF  THE  METHOD 


NUMERICAL  INTEGRATION 

Equations  (13)  are  integrated  numerically,  as  an  initial  value 
problem,  via  the  code  DVDQ  of  Krogh^^.  DVDQ  uses  a  variable  time  step, 
variable  order,  Adams  predictor-corrector  algorithm.  The  code  has  the 
advantage,  which  Is  very  Important  for  this  application  (see  p.  23),  that 
it  minimizes  the  number  of  function  evaluations  (i.e.,  air  velocity  cal¬ 
culations^®^). 

Small  Particles 

A  sufficiently  small  particle  should  essentially  follow  a 
streamline  since  both  gravity  and  inertia  effects  are  negligible.  Thus 
to  check  the  numerical  Integration  accuracy  we  computed  trajectories  of 
1  pm  diameter  water  drops  In  axi symmetric  flow  about  a  prolate  ellipsoid 
of  fineness  ratio  2  (see  Fig.  3).  Streamfunctlon  values  were  compared 
at  the  initial  (upstream)  particle  locations  and  at  the  points  where  the 
trajectories  crossed  the  extended  minor  axis  of  the  ellipsoid.  By  use 
of  the  streamf unction  gradient  along  the  extended  minor  axis,  we  deter¬ 
mined  that  the  maximum  discrepancy  in  the  axis  intersection  points  is 
0.006%.  This  is  for  the  trajectory  closest  to  the  ellipsoid  surface. 

(On  the  scale  of  the  C130A  this  represents  a  trajectory  error  of 
0.28  mm.)  Thus,  the  1  pm  water  drops  do  essentially  follow  the  stream¬ 
lines.  In  fact.  Fig.  3  is  actually  a  plot  of  the  1  pm  drop  trajectories 
used  for  this  analysis. 

Large  Particles 

Unfortunately  there  are  no  theoretical  means  nor  adequate  data 
to  check  our  calculations  for  large  particles.  The  best  that  can  be 
done  is  compare  our  calculations  with  those  Df  others.  Of  the 
results  available  in  the  literature  we  choosa  those  of  Dorsch,  Brun,  and 
Gregg^12)  as  being  the  most  useful  for  our  purposes.  They  report 


detailed  results  of  calculations  of  Impaction  of  water  drops  on  prolate 
ellipsoids  of  fineness  ratio  5.  In  our  studies,  we  find  that  calcula¬ 
tion  of  tangent  trajectories  Is  quite  sensitive  to  many  physical  and 
numerical  parameters  employed  In  the  computations.  Thus  we  have  calcu¬ 
lated  tangent  trajectories  and  compared  them  with  values  given  In  Fig.  4 
of  NACA  TN-3099. 

Dorsch  et  al.  ignore  gravity  in  their  calculations.  With  this 
simplification  we  have  axial  symmetry,  and  Eqs.  (13)  become 


3?*  =  485^  (vfx  "  vpx}  P 


dF  =  (vfr  "  vpr}  P 


(19) 


where  x  and  r  are  the  axial  and  radial  coordinates,  P  is  the  ratio  of 
Best  to  Reynolds  numbers  (Eq.  (14)),  and  S^,  variously  called:  inertial, 
impaction,  or  Stokes  number,  is 


2  pd(6/2)2  V  Pp6ZV 

sn  =  9  —  =  w 


(20) 


Here  S„  Is  defined  to  be  consistent  with  Dorsch  et  al. 

A  tangent  trajectory  grazes  the  body  surface.  For  an  axially 
symmetric  case,  any  trajectory  with  initial  r  smaller  than  rtan,  will 
impact  on  the  body.  For  fixed  body  geometry,  rtan  is  a  function  of  the 
Stokes  number,  SN>  and  a  Reynolds  number,  R^^q, 
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(21) 


Calculations  were  done  for  two  sets  of  S^,  q.  Results  are 
shown  In  Table  10.  Considering  the  difference  In  computational  approach 
(the  NACA  people  used  a  mechanical  differential  analyzer),  and  difference 
in  drag  data,  agreement  is  good. 


TABLE  10 

COMPARISON  OF  TANGENT  TRAJECTORY  RESULTS 
WITH  THOSE  OF  DORSCH  ET  AL.^12^ 


rtan 


SN 

rn,o 

Dorsch  et  al . '  ^ 

Norment  and  Zalosh 

1 

4096 

.076  i  .0014 

.0750  ±  .0001 

1/30 

512 

.0200  ±  .003 

.0150  ±  .0001 

POTENTIAL  FLOW  CALCULATIONS 

It  is  well  established  that  potential  flow  calculations 
provide  excellent  approximations  to  airflow  around  smooth  bodies  at  sub¬ 
sonic  speeds,  provided  that  the  skin  friction  boundary  layer  and  regions 
of  separated  flow  are  avoided.  Whitten computed  the  boundary  layer 
thickness  on  the  C130A  fuselage;  as  far  aft  as  FS  350",  he  found  that  the 
boundary  layer  is  no  thicker  than  3*s  inches.  Since  our  target  points  are 
greater  than  9"  away  from  the  fuselage,  boundary  layer  effects  do  not 
appear  to  be  significant. 

Hess  and  Smith  provide  ample  proof,  in  terms  of  point-by¬ 
point  comparisons  of  results,  that  their  method  for  arbitrary  three- 

H9  201 

dimensional  bodies  is  very  accurate'  *  '.We  have  compared  complete 


37.  R.  P.  Whitten,  "An  Investigation  of  Some  Aerodynamic  Factors  Affect¬ 
ing  Meteorological  Instrument  Readings  on  a  C130A  Research  Aircraft," 
Allied  Research  Associates,  AFRD  TN-6C-454  (15  May  1960). 
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trajectories  for  particles  in  analytical  potential  flow  about  an  ellip¬ 
soid,  with  potential  flow  about  a  Hess-Smlth  approximate  ellipsoid  con¬ 
sisting  of  1800  quadralateral  panels.  A  fineness  ratio  of  2,  which 
provides  a  reasonable  likeness  to  the  C130A  fuselage,  was  used.  Air 
speed,  scale,  and  air  properties  were  as  given  in  Table  1  for  the  C130A 
at  5  kft  altitude.  Results  are  shown  in  Fig.  33  for  comparison  of  tra¬ 
jectory  Intersections  with  the  extended  minor  axis.  All  of  the  Hess-Smlth 
calculation  points  are  slightly  greater  than  the  analytical  points.  The 
discrepancies  are  of  acceptable  magnitude.  The  largest  discrepancy,  for 
100  ym  drops  at  31  cm,  is  very  atypical  in  that  this  point  is  on  the  edge 
of  a  shadow  zone,  where:  trajectory  distortions  are  near  their  maxima, 
concentration  factors  become  very  large,  and  we  expect  and  find  "patho¬ 
logical"  computational  results  caused  by  trajectories  crossing  each 
other  in  this  region  of  extremely  high  concentration  gradients. 

AERODYNAMIC  PARTICLE  DRAG 

Drag  data  for  terminal  settling  of  particles  are  known  with 
high  accuracy.  Here  the  question  of  accuracy  arises  from  use  of  these 
data  for  accelerated  motion.  No  rigorous  theory  exists  for  accelerated 
motion  of  particles  at  intermediate  and  large  Reynolds  numbers  ( i . e. , 
beyond  the  Stokes  range).  In  spite  of  this,  calculations  of  the  kind 
required  here  are  routinely  done.  Fuchs  discusses  this  situation  in  his 
book^.  He  reasons  that  as  long  as  the  Reynolds  number  does  not  exceed 
several  hundred,  acceptable  results  should  be  obtained  by  use  of  steady- 
state  drag  coefficients. 

It  is  apparent  that  for  non-steady  motion  the  drag  coeffi¬ 
cient  should  be  a  function  of  acceleration  as  well  as  Reynolds  number. 
Unfortunately,  there  is  considerable  confusion  in  the  available  experi¬ 
mental  data.  Fuchs^  is  of  the  opinion  that  the  drag  should  increase 
relative  to  the  steady  state  for  accelerating  motion,  but  decrease  for 
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TRAJECTORY  ERROR  (cm) 


FIGURE  33.  Comparison  of  water  drop  trajectories  about  an  ellipsoid 
of  fineness  ratio  2  using  exact  and  approximate  potential 
airflow.  (See  text.) 
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decelerating  motion.  However,  Ingebo^33'  finds  a  decreased  drag  for 
accelerating  motion.  Ogden  and  Jayaweera^3^  recently  found  decreased 
drag  for  decelerating  motion,  though  less  in  magnitude  than  found  by 
Ingebo. 

Ingebo's  data  for  spheres  very  nicely  fit  a  power  law  relation 
between  Cp  and  R^.  His  relation  is 

CD  =  27  D'0-84  .  (22) 

Notice  the  independence  of  acceleration.  Equation  (22)  can  be  compared 
with  Stokes  law,  which  It  should  approach  In  the  limit  as  -*>  0, 

CD  =  24  R"1  ,  (23) 

and  to  an  approximation  equation  for  terminal  settling  of  spheres  at 
large  Reynolds  numbers, 

CD  =  28  R"0-85  +  0.48  .  (24) 

We  have  performed  some  exploratory  calculations  to  assess  the 
effect  of  non-steady  motion  on  concentration  factor  calculations.  These 
calculations  consist  of  two  kinds.  First,  we  subject  initially  steady¬ 
settling  water  drops  to  an  impulsive  onset  of  air  flow,  and  compute  the 
drop  response  using  (a)  steady-state  drag  coefficients,  and  (b)  Ingebo 

38.  R.  D.  Ingebo,  "Drag  Coefficients  for  Droplets  and  Solid  Spheres  in 
Clouds  Accelerating  in  Airstreams,"  NACA-TN  3762  (Sept.  1956). 

39.  T.  L.  Ogden  and  K.  0.  L.  F.  Jayaweera,  "Drag  Coefficients  of  Water 
Drops  Decelerating  in  Air,"  Quart.  J.  Roy.  Met.  Soc.  97_,  571  (1971). 


drag  coefficients*.  Results,  In  terms  of  relaxation  time**,  are  shown 
in  Fig.  34  for  the  horizontal  velocity  component  after  Impulsive  onset 
of  a  1  m/sec  air  velocity  In  the  horizontal  direction.  There  is  little 
difference  for  the  small  drops,  but  progressively  greater  dif Terence  as 
the  drop  size  increases.  Air  temperature  and  density  are  as  given  in 
Table  1  for  5  kft  altitude.  Second,  we  computed  concentration  factors 
for  a  few  representative  cases.  Results  are  shown  in  Table  11.  These 
calculations  indicate  the  possibility  of  a  significant  effect.  In  any 
case,  this  problem  needs  additional  study  because  of  the  lack  of  a  sound 
basis  for  use  of  the  steady-state  drag  coefficients  for  accelerative 
motion. 


TABLE  11 

COMPARISON  OF  STEADY-STATE  AND  INGEBO  DRAG  CONCENTRATION  FACTORS 
AT  THE  PARTICLE  REPLICATOR  SLIT  ON  THE  LOCKHEED  C130A 

(Conditions  are  as  given  In  Table  1  for  5  kft  altitude) 


Water  drop 
diameter 
(Pm) 

Concentration 

Factors 

Steady  State 

Ingebo 

100 

1.40 

1.24 

150 

1.33 

1.10 

*  Numerical  problems  arise  wh?n  Ingebo's  relation  is  used  at  very  low 
Reynolds  numbers.  Therefor?,  we  switch  over  to  Stokes'  relation 
where  the  Stokes  and  Ingebo  curves  cross,  at  =  0.47896. 


**  Relaxation  time  is  time  reqjired  for  a  particle  to  reach  1-1/e  of 
its  final  velocity. 
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RELAXATION  TIME  (sec) 


WATER  DROP  DIAMETER  (ym) 


FIGURE  34.  Relaxation  times  of  water  drops,  computed  by  steady-state  and 
Ingebo  drag  coefficients,  in  response  to  a  1  m/sec  impulsive 
horizontal  flow. 


SENSITIVITY  TO  FLUX  TUBE  STRUCTURE 

It  might  be  supposed  that  concentration  factor  accuracy  is 
sensitive  to  the  structure  of  the  particle  flux  tube,  since  the  flux 
tube  is  used  to  approximate  an  Infinitesimal  property  of  the  flow. 
Accordingly,  sensitivity  studies  were  done  for  water  drops  at  the  C130A 
particle  replicator  at  5  kft  altitude.  Results  are  shown  in  Table  12. 
Over  a  wi^  .■  of  tube  structures,  maximum  differences  of  less  than 
2%  are  fou.id.  These  differences  are  no  greater  than  those  found  with 
the  same  tube  structure,  but  with  different  initialization  guesses. 

PARTICLE  SHAPE  AND  ORIENTATION 

We  assume  spherical  shape  for  water  drops.  This  simplification 

is  consistent  with  the  observation  of  negligible  shape  deformation  for 

terminally  settling  drops  smaller  than  ^400  ym  diameter.  Furthermore, 

the  smaller  drops  are  not  expected  to  deform  appreciably  in  the  weak 

shear  found  outside  the  boundary  layer.  Since  concentration  factors 

approach  unity  for  drops  larger  than  400  ym,  we  conclude  that  for  water 

drops  we  have  neither  a  shape  nor  orientation  problem. 

For  ice  columns  we  have  computed  concentration  factors  for 

"axis-on"  and  "end-on"  column  orientations  relative  to  the  drag  vector 

(see  Fig.  10).  Many  have  observed  that  the  "axis-on"  orientation  is 

(341 

preferred,  and  Bragg  et  al.  '  observe  that  this  orientation  is  very 
rapidly  attained.  Thus,  for  well  formed,  unrimed  ice  columns,  our 
"axis-on"  concentration  factor  results  can  be  accepted  with  confidence. 

On  the  other  hand,  fragmented  or  rimed  columns  tend  to  oscillate, 
tumble,  or,  In  extreme  cases,  settle  with  their  long  axes  oriented 
vertically.  Our  results  may  not  be  representative  for  such  cases. 

As  yet  we  have  not  applied  our  method  to  other  ice  forms: 
plates,  aggregates,  dendrites,  etc.  It  will  be  interesting  to  see  how 
results  for  these  forms  compare  with  results  for  drops  and  columns. 
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TABLE  12 


SENSITIVITY  OF  Cf  TO  FLUX  TUBE  STRUCTURE 

100  pm  Water  Drops  at  the  C130A  Replicator  Slit 
5  kft  altitude 


A.  Number  of  Trajectories  on  the  Tube  Surface  (r  =  .01,  e  =  .2) 

n 

N  CF 

4  1.399 

6  1.388 

8  1.391 

B.  Target  Window  Radius  and  Convergence  Tolerance 


Cf 

e 

rw  -  .005 

.01 

.02 

0.05 

0.1 

1.390 

1.394 

1.413 

0.2 

0.4 

1.393 

1.400 

1.412 
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SUMMARY  AND  DISCUSSION 


We  are  satisfied  that  integration  and  flow  calculation  numeri¬ 
cal  errors  are  negligible.  The  digital  descriptions  of  the  Lockheed 
airplane  fuselages  are  very  accurate  (^.02  Inch).  For  the  Cessna,  some 
interpolation  was  necessary,  but  Individual  coordinates  are  known  to 
better  than  2*Hnch  accuracy.  Provided  that  the  sKin  friction  boundary 
layer  and  regions  of  separated  flow  are  avoided,  the  potential  flow 
calculations  yield  accurate  simulation  of  the  air  flow.  Steady-state 
drag  for  spheres  and  columns  is  accurately  known,  well  beyond  the  range 
of  Reynolds  numbers  encountered  (Appendix  B). 

For  water  drops,  the  only  major  source  of  error  appears  to  be 
use  of  the  steady  drag  data  for  accelerative  motion.  Sensitivity 
results  presented  above  (i.e.,  comparisons  of  steady  drag  results  with 
results  using  Ingebo's  drag  data)  would  appear  to  overstate  the  error  in 
light  of  the  results  of  Ogden  and  Jayaweera^).  We  conclude  that  con¬ 
centration  factor  errors  for  water  drops  are  less  than  20%,  except  near 
the  edge  of  a  shadow  zone.  We  suggest  use  of  the  method  presented  here 
to  correct  existing  data. 

For  ice,  the  situation  is  much  more  complex.  On  the  basis  of 
the  error  analyses  above,  and  the  results  presented  in  the  preceding 
chapter,  we  might  claim  an  accuracy  of  30-40%  for  columns.  On  the  other 
hand,  there  are  a  number  of  uncertainties  that,  are  essentially  beyond 
our  control.  For  example,  our  results  strictly  apply  only  to  well 
formed,  unrimed  ice  crystals.  We  have  used  typical  column  dimensions 
and  densities,  but  these  properties  of  ice  crystals  are  observed  to  vary 
widely.  For  the  Knollenberg  spectrometers,  the  question  of  how  to 
relate  instrument-registered  ice  particle  dimensions  with  true  particle 
shape,  dimensions,  and  density  is  an  open  one.  We  must  defer  assignment 
of  a  quantitative  error  estimate  until  other  crystal  forms  are  studied. 
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CONCL JSIONS 


A  method  has  been  developed  to  compute  concentration  distor¬ 
tion  of  particulates  of  arbitrary  size  caused  by  airflow  around  airplane 
fuselages.  The  method  combines  a  three-dimensional  potential  airflow 
solution  with  three-dimensional  hydrometeor  trajectories  around  the  air¬ 
plane.  It  has  been  applied  to  hydrometeor  sampling  instruments  on  three 
cloud  physics  research  airplanes:  Lockheed  C130A  and  C130E  transports, 
and  a  Cessna  Citation  jet.  Results  have  been  obtained  for  water  drops 
and  ice  columns. 

On  a  Lockheed  C130A  a  formvar  particle  replicator  mounted  just 
aft  of  the  cockpit  was  studied.  For  both  water  drops  and  ice  columns, 
substantial  concentration  enhan  :ement  (from  30-60%)  over  a  wide  range  of 
particle  sizes  was  found. 

On  the  Lockheed  C130E  a  formvar  particle  replicator  mounted  in 
the  forward  cargo  compartment  wis  studied.  At  5  kft  altitude,  the  calcu¬ 
lation  indicated  that  water  drops  over  a  range  from  about  90-300  pm 
diameters  cannot  be  sampled  by  the  instrument.  For  other  drop  sizes 
severe  concentration  enhancement  is  Indicated.  This  instrument  is  being 
redesigned  to  extend  the  sampling  point  8  inches  further  from  the  fuse¬ 
lage.  At  the  new  point  the  cal  :ulations  indicate  concentration  enhance¬ 
ment  (up  to  40%)  over  a  substantial  range  ol  sizes. 

Knoll enberg  particle  .pectrometer;  mounted  on  the  emeracncy 
exit  door  of  the  Cessna  Citatio  i  were  studied.  At  the  precipitation 
particle  spectrometer,  concentration  enhancement  (up  to  80%)  over  a  wide 
range  of  sizes  for  water  drops  < nd  ice  columns  are  indicated.  For  the 
cloud  particle  spectrometer,  water  drops  over  a  narrow  range  near  100  pm 
diameter  cannot  be  sampled,  nor  can  ice  greater  in  size  than  about 
300  pm  for  solid  columns  and  800  pm  for  hollow  columns.  Severe  concen¬ 
tration  enhancement  is  indicated  for  other  sizes. 

A  thorough  error  analysis  indicates  concentration  factor 
accuracy  of  *20%  for  water  drops.  We  suggest  use  of  the  method  to  cor¬ 
rect  existing  water  drop  data.  Accuracy  assessment  of  ice  particle 
results  is  deferred  until  additional  crystal  forms  are  studied. 
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APPENDIX  A.  ITERATIVE  PROCEDURES  TO  DETERMINE  THE  PARTICLE 
TRAJECTORY  THROUGH  A  POINT  IN  SPACE 


If  the  flow  field  is  so  complicated  that  particle  trajectories 
must  be  calculated  numerically,  the  particle  velocity  through  a  point 
must  be  known  to  determine  the  particle  trajectory  through  that  point. 
Since  in  this  work  we  do  not  have  a^  priori  knowledge  of  particle  velo¬ 
cities  at  target  points,  we  cannot  determine  trajectories  by  backward 
calculation  in  time.  Instead,  we  must  begin  at  a  point  far  upstream 
where  the  particle  velocity  is  known,  and  calculate  forward  in  time. 

Since  we  also  do  not  have  a  priori  knowledge  of  the  initial  point  of 
the  desired  trajectory,  the  trajectory  must  be  determined  by  a  trial  - 
and-error  iterative  process.  The  tneory  is  given  in  the  Concentration 
Factor  section  (p.  13).  Here,  we  develop  the  detailed  equations  and 
describe  their  application. 

Our  fundamental  equation  is  a  first  order  approximation  to  a 
Taylor's  series  that  relates  coordinate  points  in  the  initial  and  tar¬ 
get  planes,  i.e.,  planes  perpendicular  to  the  particle  velocity  vector 
far  upstream  and  at  the  target.  Let  our  target  point  be  Zt  Q,  and  Zt  be 
an  arbitrary  point  in  the  target  plane  which  is  close  to  Zt  Q.  Let  Z  be 
the  point  of  intersection  in  the  initial  plane  of  the  particle  trajectory 
through  Z^.  Our  fundamental  equation  is 

Z  =  B  +  C(Zt  -  Zt>0)  ,  (A-l) 

where  all  quantities  are  complex,  e.g.,  and  B  =  b^+ib^. 

(st  and  denote  the  Cartesian  coordinates  of  a  point  in  the  target 
plane). 

We  are  given  Zt  q  and  we  begin  our  iterative  procedure  by 
guessing  two  sets  of  values  for  Z^  and  Z,  which  we  label  Z^  y  Z1  and 
Zf  2»  ^2»  In  terms  of  the  real  and  imaginary  parts  of  Eq.  (A-l)  these 
points  yield  four  simultaneous  equations  in  the  four  unknowns,  b  ,  b  , 

S>  S 

c^,  c^  which  in  matrix  notation  is 
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Preceding  page  blank 


1 

0 

4tt,l 

b 

c 

C1 

0 

1 

4Et,l 

4Et,l 

• 

b£ 

*1 

1 

0 

4tt,2 

-1Et.2 

<2 

0 

1 

«t,2 

4H,2 

52 

(A-2) 


where  ActJ  =  5t>1  -  Ctf0.  A€t#1  -  Ct§1  -  «t(0.  etc.  This  system  of 

equations  is  solved  for  b  and  b?,  which,  as  seen  from  Eq.  (A-l),  pro¬ 
vide  the  next  estimate  of  the  desired  initial  coordinates  Z.  Using 
these  initial  coordinates,  which  we  label  Z3>  a  trajectory  is  computed 
and  the  corresonding  3  1 >  determined.  Then  we  repeat  the  matrix 
inversion  using 


Z2  '  B  +  CiZt,2 
Z3  3  B  +  C4Zt,3 

and  determine  a  revised  B,  and  so  on. 

This  procedure  wcrks  well  but  it  uses  only  the  most  recent 
two  sets  of  results.  If  tie  method  has  not  converged  after  calculation 
of  two  trajectories,  we  sw  tch  over  to  a  least  squares  procedure  which 
uses  all  trajectory  data  that  have  been  generated.  The  real  and  imaginary 
parts  of  Eq.  (A-l)  are 


5  3  bc  +  Cc  45t  '  Cc  4Ct 


£  =  b?  +  c5 


After  eliminating  the  term  in  c?,  we  obtain 


(A-3) 
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(A-4) 


R  = 


bt  4Ct  *  b?  act  +  ct|4Zt|2 


where 


R  =  ;  +  A£t  K 


(A-5) 


b  ,  b  ,  and  c  are  determined  such  that 

C  »  £ 

l(R  -  bt  b5u  -  cc|4Zt|2)2 

is  minimized,  where  the  summation  is  over  all  of  the  trajectory  results. 
The  normal  equations  matrix  is 


[Uct)2  l  Act  ACt  l  ACt|AZt|2 

b 

C 

I  «t  A 

K aO2  I  A5t|AZt  2 

• 

bE 

= 

l  A 

l[W2)2  _ 

cc 

IlAZtl2R 

Iteration  is  continued  until  the  current  Z. ,  labeled  Z.  . ,  is  within  a 

t  tfj 

prescribed  tolerance  (Eq.  (11))  of  being  equal  to  Zt  q. 

In  an  effort  to  improve  convergence  efficiency,  a  second- 

order  version  of  the  iteration  procedure  was  tried.  Our  fundamental 
equation  becomes 


2  -  8  +  C(Zt  -  Ztj0)  ••  0(Zt  -  Ztj0)2  .  (A-7) 

This  yielded  no  improvement  in  convergence,  and  was  abandoned  in  favor 
of  the  simpler  first  order  method. 
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A  code  has  been  developed  that  performs  automatically  the 
calculations  required  to  determine  a  concentration  factor.  The  code 
Is  supplied  with  flight  conditions  and  the  output  of  the  Hess-Smith 
code  so  that  velocities  at  arbitrary  points  in  space  can  be  computed. 

In  addition,  it  Is  given  the  ta 'get  point  coordinates,  the  target  plane 
window  radius,  r  ,  the  tolerance,  e,  the  number  of  trajectories  to  be 

W 

used  to  define  the  flux  tube,  and  two  sets  of  initial  and  target  plane 
coordinates.  The  code  first  finds  the  trajectory  through  the  central 
target  point  by  iteration  until  Eq.  (11)  is  satisfied.  (In  the  event 
that  convergence  is  not  achieved  in  25  attempts  the  concentration  factor 
calculation  is  abandoned.)  Then  it  finds  each  of  the  flux  tube  trajec¬ 
tories  in  a  similar  manner.  It  computes  the  areas  of  the  polygons 
defined  by  the  trajectory  intersections  in  the  initial  and  target  planes, 
and  computes  the  concentration  factor  as  the  ratio  of  the  areas. 

In  the  course  of  the  iterative  calculations,  impaction  on  the 
airplane  body  frequently  occurs.  In  this  event  the  initial  coordinates, 
c,  5,  are  incremented  by  1%  of  their  values  such  as  to  shift  the  tra¬ 
jectory  away  from  the  body,  and  the  calculation  is  begun  again.  If  a 
total  of  25  impactions  are  registered  in  attempting  to  compute  the 
trajectory  through  a  target  point,  the  concentration  factor  calculation 
is  abandoned. 
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APPENDIX  B.  AERODYNAMIC  DRAG  ON  MOVING  PARTICLES 


GENERAL  CONSIDERATIONS 

In  this  work,  we  take  particle  motion  relative  to  the  air 
to  result  from  two  phenomena:  (1)  free-fall,  gravity  settling,  and 
(2)  particle  Inertia  relative  to  air  being  displaced  by  a  passing  air¬ 
plane.  Turbulence  effects  are  unimportant  over  the  time  intervals 
considered  here. 

An  adequate  data  base  exists  to  allow  accurate  calculation  of 
terminal  settling  of  particles  of  simple  shapes.  For  accelerative 
motion  of  particles,  there  are  neither  adequate  data  nor  theory.  Never¬ 
theless,  we  use  the  terminal  settling  drag  data  for  accelerative  motion 
computations.  This  situation  is  discussed  in  the  chapter  on  Accuracy. 

It  is  easy  to  show  from  Eq.  (12)  that  the  Best  number  for 

p 

terminal  particle  settling,  BN  j  (=  R^  j  Cp) ,  is 

Bnt-^  .  (B-l) 

nZAp 

Symbols  are  as  defined  previously  (also  see  Appendix  D).  Notice  that 
j  is  independent  of  the  particle  settling  speed.  Since  T  is  a 
function  of  all  of  the  particle  and  fluid  properties  that  determine 
flow  around  the  particle,  the  settling  speed  must  be  a  function  of  B^  j. 
Thus 

RN,T  =  fl(BN,T}  ’  (B"2) 

where  R^  y,  the  Reynolds  number  for  terminal  settling  is 
p  Vt6 

n  _  1  fo  o\ 
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Data  needed  to  determine  f^(BN  y)  are  available  in  the  literature,  and 
are  discussed  below. 

In  computing  particle  trajectories  through  complex  flow  fields, 
we  integrate  Eqs.  (13)  numerically,  step-by-step  forward  in  time.  At 
any  time,  and  corresponding  point  in  space,  the  particle  Reynolds  number 
is  known.  We  must  determine  P  -  B^/R^  from  the  Reynolds  number.  For 
this  purpose,  we  use  the  same  data  as  Is  used  to  determine  function 
fl(BN  y),  but  in  reverse  order,  to  determine  the  function,  f^ 

Bn  =  f2(RN)  •  (B-4) 

As  noted  previously,  f£  should  be  a  function  of  acceleration  as  well  as 
Reynolds  number,  but  for  moderate  Reynolds  numbers  and  accelerations 
Eq.  (B-4)  seems  to  yield  an  ac equate  approximation. 

In  the  remainder  of  this  Appendix,  we  discuss  our  determina¬ 
tions  of  the  functions  f1  and  f 


SPHERES 


The  Best  number  for  terminal  settling  of  spheres  reduces  to 


4 

3 


3 

ppp6  g 

“V" 


(B-5) 


where  pp  is  sphere  density, 
a  large  body  of  experimental 


Davies fits  the  following  functions  to 
data: 


40.  C.  N.  Davies,  "Definitive  Equations  for  the  Fluid  Resistance  of 
Spheres,"  Proc.  Phys.  Soc  (London)  57,  259  (1945). 
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rnj  »  -!jf  -  2. 3363x1 0"4  B2>T  +  2. 01 54x1 0"6  B3  T 

-  6 .9105x1 0”^  bJ>t  ,  BN  T  <  140  (B-6) 

l°g10  rn  t  =  “  1,29536  +  0.986  (log]0  B^-j.) 

-  0.046677  (log1Q  BNJ)2 
+  0.0011235  (log1Q  BNjT)3  , 

100  <  BN  t  <  4.5xl07  (B-7) 

To  find  f2,  we  have  fit  the  data  tabulated  by  Davies  to 
polynomials  in  RN  j,  viz. 

n 

vl  vi  *  ^B_8^ 

J-o 

The  polynomial  coefficients  determined  by  least  squares  are  given  in 
Table  13. 
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TABLE  13 


POLYNOMIAL  COEFFICIENTS  RELATING  BEST  NUMBER 
TO  REYNOLDS  NUMBER  FOR  SPHERES 


BN 


n 


I  ‘A 


Reynolds  No.  Range 
0.05  <  <  3 


0 

1 

2 

3 


a . 

_ J_ 

0 

24.167 

3.2540 

-0.23564 


3  <  RN  <  330 


0  -28.339 

1  38.969 

2  0.73204 

3  -0.00056084 


330  < 


0  0 

1  93.462 

2  0.37576 
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COLUMNS 


Experimental  data  for  "axison",  terminal  settling  of  cylinders 
is  reported  by  Kajikawa^4^,  Jayaweera  and  Mason^33^,  and  Jayaweera  and 
Cottis^42^.  Though  these  data  are  for  circular  based  cylinders,  Jayaweera 
and  Cottis^42^  maintain  ihat  they  can  be  used  for  hexagonal  based  cylin¬ 
ders  without  significant  error.  For  terminal  setting,  the  Best  number  is 


(B-9) 


where  5  and  i  are  as  defined  in  Fig.  10,  and  the  Reynolds  number  is 


fB-10) 


By  use  of  a  Kargel  Reflecting  Projector,  we  transferred  the 
graphically  displayed  data  point?  of  Kajikawa  and  Jawaweera  and  Cottis 
to  lined  graph  paper.  Each  datun  was  independently  recorded  at  least 
twice.  With  two  exceptions  thes?  authors  made  their  measurements  for 
the  same  6/t  values.  The  data  set  for  one  exceptional  6/1  (0.2)  is 
inconsistent  with  the  remainder  and  was  discarded.  Data  for  6/i=0, 
measured  by  Jayaweera  and  Mason^33^,  were  taken  from  the  presentation 
of  Jayaweera  and  Cottis.  The  d;  ta  from  the  two  papers  were  merged  and 
fitted  via  least  squares  to  polynomial  functions  of  log10  (R^  T)  vs. 
logiO  (Bn  j)-  These  polynomials  are  described  in  Table  14.  In 
Table  15  are  the  least  squares  coefficients  for  the  inverse  polynomials 
used  to  evaluate  log-j q  B^  from  log-jQ  R^. 


41.  M.  Kajikawa,  "A  Model  Experimental  Study  on  the  Falling  Velocity  of 
Ice  Crystals,"  J.  Met.  Soc.  Japan  4£,  367  (1971). 

42.  K.  0.  L.  F.  Jayaweera  and  R.  E.  Cottis.  "Fall  Velocities  of  Plate- 

Like  and  Columnar  Ice  Crystals,"  Quart.  J.  Roy.  Met.  Soc.  95,  703 
(1969).  ~ 
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TABLE  14 


POLYNOMIAL  COEFFICIENTS  RELATING  REYNOLDS  NUMBER  TO  BEST  NUMBER 
FOR  AXIS-ON  TERMINAL  SETTLING  OF  CYLINDERS 


4 

,o910RN.T  ■  I  bj(1O910BN,T)J 
j*0 


Number  of 
data  points 


6/1 

V 

bi 

b2 

b3 

b4 

J&C 

Kaj 

Total 

0.0 

-  .76690 

.85642 

-.060794 

.011948 

-.0019287 

26 

0 

26 

0.1 

-  .84486 

.87235 

-.065890 

.013440 

-.0018927 

20 

19 

39 

0.5 

-1.1073 

.98200 

-.054802 

0 

0 

17 

16 

33 

1.0 

-1.3140 

1.0115 

-.035804 

-.0034749 

0 

21 

16 

37 

6/£ 

Range  n  R^  j 

Range  in  j 

0.0 

.188  -  84 

1.11  -  8400 

0.1 

.0219  -  71 

.136  -  630 

0.5 

.00664  -  107 

.090  -  12600 

1.0 

.00398  -  98 

.088  -  12300 
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TABLE  15 


POLYNOMIAL  COEFFICIENTS  RELATING  BEST  NUMBER  TO  REYNOLDS  NUMBER 

FOR  AXIS-ON  CYLINDER  DRAG 


1og10BN  ■  I  aj  [,°S10Rn] 
j*0 

3 

1 Zl 

0.0 

a0 

31 

a2 

a3 

a4 

.95007 

1.3137 

.073827 

-.023087 

.024530 

0.1 

1 .0320 

1.2841 

.090954 

.0064871 

0 

0.5 

1 .2048 

1.1746 

.10077 

.015304 

0 

1.0 

1.3679 

1.1124 

.087279 

.018643 

0 
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Kajlkawa^41)  gives  data  for  "end-on"  terminal  settling  of 
columns.  As  noted  previously,  we  have  good  reason  to  reject  this 
settling  orientation  for  ice  columns  In  the  free  air.  On  the  other 
hand,  there  Is  some  probability,  albeit  small,  that  columns  will  orient 
with  their  axes  parallel  with  the  drag  vector.  Therefore,  we  have 
reduced  Kajlkawa's  end-on  drag  data  as  described  above  to  yield  the 
polynomials  given  In  Table  16.  Concentration  factor  results  shown  in 
Figs.  15-18  and  tabulated  In  Table  6  for  end-on  column  orientation  were 
computed  using  these  polynomials.  They  were  used  to  compute  drag  on 
the  columns  during  the  trajectory  calculations,  though  settling  was 
computed  via  use  of  the  axis-on  polynomials. 


TABLE  16 

POLYNOMIAL  COEFF  CIENTS  RELATING  BEST  NUMBER  TO 
REYNOLDS  NUMBLR  FOR  END-ON  CYLINDER  DRAG 

,09lOBN  ‘  I  8J<'°9lORN)J 
j*0 


6/ 1 

aO 

al 

a2 

Number  of 
data  points 

1/3 

1.6098 

1 .0259 

0.0 

3 

1/2 

1.4616 

0.99954 

0.0 

3 

1/1 

1.4513 

1 .0724 

0.021773 

4 
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APPLICATION  OF  THE  DRAG  POLYNOMIALS 

Given  drop  diameter  or  crystal  dimensions,  particle  density 
and  atmospheric  properties  ,  y  Is  calculated  by  Eq.  (B-5)  or  (B-9). 
For  watar  drops,  R^  y  is  calculated  from  Eq.  (B-6)  or  (B-7).  For  Ice 
columns,  log1Q  Rfj  j  values  are  calculated  via  the  two  polynomials  In 
Table  14  that  bracket  the  given  5/1  value.  Log.|Q(RN  y)  for  the  column 
Is  found  by  linear  Interpolation  In  6/1.  Finally,  Vy  Is  computed  as 


V 


T 


nRN,T 

p6 


(B-ll ) 


To  compute  the  drag  on  a  particle  during  Integration  of  the 
particle  equations  of  motion,  a  similar  procedure  Is  followed.  The 
Reynolds  number  Is  computed  from  Eq.  (16).  For  a  water  drop,  BN  Is 
computed  from  Eq.  (B-8).  For  Ice  columns,  1  og-j q  Bn  Is  computed  via 
the  polynomials  In  Table  15  or  16  that  bracket  the  given  6/fc  value. 
Finally,  log1Q  BN  for  the  given  6/t  Is  obtained  by  linear  interpolation 
In  <5/1. 
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APPENDIX  C.  PROJECTED  DIMENSION  0:  RANDOMLY  ORIENTED  RECTANGLES 


Here  we  wish  to  relate  actual  Ice  column  dimensions,  6  and 
to  the  mean  characteristic  dimension  recorded  by  the  Knollenberg  specto- 
meters.  We  assume  that  Ice  columns  settle  with  their  long  axes  at  ran¬ 
dom  attitudes  In  the  horizontal  plane.  A  horizontally  oriented  linear 
sizing  instrument  will  register  a  projection  of  the  column  dimensions 
onto  its  linear  axis. 

A  plan  view  of  the  geometry  Is  shown  In  Fig.  C.l,  where  the 
long  axis  of  length  *  makes  an  angle  e  with  the  projection  direction. 

The  short  axis  has  length  6.  From  the  drawng  we  see  that  the  projected 
dimension,  v.  Is 

v  =  i  sin  e  +  6  cos  e  .  (C-l) 

The  probability  that  e  will  assume  a  value  between  e  and 
e  +  de  Is  simply 


dp(e)  *  ^  de  ,  0  <  e  <  | 


(C-2) 


Therefore,  for  fixed  i  and  6  the  mean  value  of  v  is 


<v> 


TT  /  2  w/2 

/  7  *  f 1  sin  e  +  6  cos  e)  de 


<v>  ■-(*  +  «) 

n 


(C-3) 
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Preceding  page  blank 


APPENDIX  D.  GLOSSARY  OF  SYMBOLS 

cross-sectional  area  of  a  particle  flux  tube  in  the  free- 
o 

stream  (m  ) 

p 

area  of  particle  projected  in  its  direction  of  motion  (m  ) 

cross-sectional  area  of  a  particle  flux  tube  in  the  target 
plane  (nr) 

Best  number 

coordinate  perpendicular  to  fuselage  symmetry  axis  (English 
units) 

drag  coefficient  (dimensionless) 
concentration  factor  (Eq.  (1)) 
particle  concentration  ratio  (Eq.  (2)) 
particle  flux  In  the  free-stream  (kg/(m  -sec)) 

Froude  number 

2 

particle  flux  at  a  target  point  (kg/(m  -sec)) 

coordinate  parallel  with  fuselage  axis  (English  units) 

2 

gravity  acceleration  constant  (9.8  m/sec  ) 

2 

acceleration  of  gravity  (m/sec  ) 

/T 

ice  column  length  (ym) 

characteristic  dimension  of  an  airplane  (m,  feet  or  inches) 
particle  mass  (kg) 

mass  transfer  rate  through  a  particle  flux  tube  (kg/sec) 
ratio  of  Best  to  Reynolds  numbers  (P  =  B^/R^) 
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p 


T 


r 


w 


fx’  vfy’  vfz 


V*  Vpy’  V 


V 

-► 

V, 


T 

ML 


x,  y,  z 


Z 


D 


P 


P 


ratio  of  Best  to  Reynolds  numbers  for  terminal  particle 
settling 

radius  of  particle  flux  tube  cross  section  In  the  target 
plane  (dimensionless) 

Reynolds  number 

Stokes  number 

time  (sec) 

air  velocity  components  (dimensionless) 

particle  velocity  components  (dimensionless) 

air  speed  at  the  target  point  (dimensionless) 

terminal  settling  speed  of  a  particle  (dimensionless) 

free-stream  air  speed  (m/sec) 

air  velocity  (m/sec) 

particle  velocity  (m/sec) 

air  speed  at  the  target  point  (m/sec) 

terminal  settling  speed  of  a  particle  (m/sec) 

vertical  coordinate  in  fuselage  system  (English  units) 

space  coordinates  (dimensionless) 

complex  point  in  a  plane  perpendicular  to  a  particle 
flux  tube  (dimensionless)  (Z  =  c  +  U) 

target  point  (Zt0  =  ttj0  +  Ut>0) 

imaginary  coordinate  of  a  complex  point  in  a  plane  per¬ 
pendicular  to  a  particle  flux  tube  (dimensionless) 

3 

air  density  (kg/m  ) 
particle  density  (kg/m  ) 
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time  (dimensionless) 

particle  dimension  ( :-n) .  (water  drop  diameter  or  column 
base  width) 

fraction  of  rw  used  for  trajectory  conversion  criterion 
(see  Eq.  (11)) 

real  coordinate  of  a  complex  point  In  a  plane  perpendlcu 
lar  to  a  particle  flux  tube  (dimensionless) 

air  viscosity  (kg/(m-sec)) 

projected  dimension  of  an  Ice  column  (pm) 


-  97  - 


REFERENCES 


1.  P.  Spyers-Duran  and  R.  R.  Braham,  "An  Airborne  Continuous  Cioud 
Particle  Replicator,"  J.  Appl .  Meteor.  6,  1108  (1967). 

2.  J.  Hallett,  R.  W.  Hanaway,  and  P  B.  Wagner,  "Design  and  Construction 
of  a  New  Cloud  Particle  Repll cater  for  Use  on  a  Pressurized  Aircraft," 
Desert  Research  Institute,  Reno,  Nevada,  AFCRL-72-0410  (31  May  197?). 
AD- 753  091 . 

3.  R.  G.  Knoll enberg,  "The  Optical  Array:  An  Alternative  to  Scattering 
or  Extinction  for  Airborne  Particle  Size  Determination, "  J.  Appl. 
Meteor.  9,  86  (1970). 

4.  R.  G.  Dorsch  and  R.  J.  Brun,  "Variation  of  Local  Liquid-Water  Concen¬ 
tration  About  an  Ellipsoid  of  Fineness  Ratio  5  Moving  In  a  Droplet 
Field,"  NACA-TN-31 53  (July  1954). 

5.  G.  F.  Carrier,  M.  Krook,  C.  E.  Pearson,  Functions  of  a  Complex 
Variable  (McGraw-Hill  Book  Company,  1966]"!  Sec.  2  6. 

6.  N.  A.  Fuchs,  The  Mechanics  of  Aerosols,  Translated  by  R.  E.  Daisley 
and  M.  Fuchs,  edited  by  t.  N.  bavles  (MacMillan,  New  York,  1964). 
Chapter  III. 

7.  C.  N.  Davies  and  M.  Aylward,  "The  Trajectories  of  Heavy,  Solid 
Particles  In  a  Two-Dimensional  Jet  of  Ideal  Fluid  Impinging  Normally 
upon  a  Plate,"  Proc.  Phys.  Soc.  (London)  B64,  889  (1951). 

8.  C.  N.  Davies  and  C.  V.  Peetz,  "Impingement  of  Particles  on  a  Trans¬ 
verse  Cylinder,"  Proc.  Roy.  Soc.  (London)  A234,  269  (1956). 

9.  K.  V.  Beard  and  S.  N.  3rover,  "Numerical  Collision  Efficiencies  for 
Small  Raindrops  Colllcing  with  Micron  Size  Particles,"  J.  Atmos. 

Sci.  3 1,  543  (1974). 

10.  R.  L.  Pltter  and  H.  R.  Pruppacher,  "A  Numerical  Investigation  of 
Collision  Efficiencies  of  Simple  Ice  Plates  Colliding  with  Super¬ 
cooled  Water  Drops,"  J.  Atmi's.  Sci.  31_,  551  (1974). 

11.  I.  Langmuir  and  K.  B.  Blodgett,  "Mathematical  Investigation  of  Water 
Droplet  Trajectories,"  General  Electric  Company,  Report  RL-225  (1945). 

12.  R.  G.  Dorsch,  R.  J.  Brun,  and  J.  L.  Gregg,  "Impingement  of  Water 
Droplets  on  an  Ellipse  id  with  Fineness  Ratio  5  in  Axisymmetric  Flow," 
NACA  TN-3099  (March  1954).  (Also  see  NACA  TN-2952,  2999,  3047,  3155, 
3410,  3153,  3586). 


-  99  - 


Preceding  pege  blank 


13.  H.  G.  Norment,  “Research  on  Circulation  In  Nuclear  Clouds,  II,"  Technical 
Operations,  Inc.,  Report  TO-B  64-102  (November  1964).  AD-361  074. 

14.  B.  Etkln,  "Interaction  of  Precipitation  with  Complex  Flows,"  Proceed¬ 
ings  of  the  Third  Intei  national  Conference  on  Wind  Effects  on 
Buildings  and  Structures,  (tokyo,  1971). 

15.  S.  A.  Morsi  and  A.  J.  Alexander,  "An  Investigation  of  Particle  Tra¬ 
jectories  In  Two-Phase  Flow  Systems,"  J.  Fluid  Mech.  55,  193  (1972). 

16.  S.  Kuo-Kai  Chal ,  "Dropiet  Trajectories  Around  Aircraft  Wing,"  Thesis, 

U.  Nevada,  Reno  (November  1973). 

17.  F.  T.  Krogh,  "Variable  Order  Integrators  for  Numerical  Solution  of 
Ordinary  Differential  :quations,"  Jet  Propulsion  Lab  Technology 
Utilization  Document  N).  CP-2308  (November  1970). 

18.  T.  E.  Hull,  W.  H.  Enrl  |ht,  B.  M.  Fellen,  and  A.  E.  Sedgwick,  "Com¬ 
paring  Numerical  Metho  is  for  Ordinary  Differential  Equations," 

SIAM  J.  Numer.  Anal.  9.  603  (1972). 

19.  J.  L.  Hess  and  A.  M.  0  Smith,  "Calculation  of  Non-Lifting  Potential 
Flow  about  Arbitrary  Ti  ree-DImenslonal  Bodies,"  McDonnell  Douglas 
Report  E.  S.  40622  (15  March  1962).  AD-282  255. 

20.  J.  L.  Hess  and  A.  M.  0  Smith,  "Calculation  of  Potential  Flow  About 
Arbitrary  Bodies,"  in  f rogress  in  Aeronautical  Sciences,  Vol .  8, 
edited  by  D.  Kuchemann  (Pergammon  Press,  New  York,  1967). 

21.  F.  A.  Woodward,  "Analysis  and  Design  of  Wing-Body  Combinations  at 
Subsonic  and  Supersonic  Speeds,"  J.  Aircraft  5,  528  (1968). 

22.  R.  Gunn  and  G.  0.  Kinzer,  "The  Terminal  Velocity  of  Fall  for  Water 
Droplets  in  Stagnant  Air,"  J.  Meteor.  6,  243  (1949). 

23.  K.  V.  Beard  and  H.  R.  Pruppacher,  "A  Determination  of  the  Terminal 
Velocity  and  Drag  of  Small  Water  Drops  by  Means  of  a  Wind  Tunnel," 

J.  Atmos.  Sci.  26,  1066  (1969). 

24.  H.  R.  Pruppacher  and  K.  V.  Beard,  "A  Wind  Tunnel  Investigation  of 
the  Internal  Circulation  and  Shape  of  Water  Drops  Falling  at  Terminal 
Velocity  in  Air,"  Quart.  J.  Roy.  Met.  Soc.  96,  247  (1970). 

25.  H.  R.  Pruppacher  and  R.  L.  Pitter,  "A  Semi -Empirical  Determination 
of  the  Shape  of  Cloud  and  Rain  Drops,"  J.  Atmos.  Sci.  28,  86  (1971). 

26.  C.  Magono  and  C.  W.  Lee,  J.  Fac.  of  Sci.,  Hokkaido  U.,  Ser.  VII, 

Vol.  II,  321  (1966). 


-  100  - 


27.  N.  H.  Fletcher,  The  Physics  of  Ra Inclouds  (Cambridge  University 
Press,  1966). 

28.  B.  J.  Nason,  The  Physics  of  Clouds  (Clarendon  Press,  Oxford,  1971). 

29.  A.  Ono,  "Growth  Node  of  Ice  Crystals  In  Natural  Clouds,"  J.  Atmos. 

Scl .  27,  649  (1970). 

30.  A.  J.  Heymsfleld  and  R.  G.  Knollenberg,  "Properties  of  Cirrus 
Generating  Cells,"  J.  Atmos.  Scl.  2£,  1358  (1972). 

31.  A.  H.  Auer  and  D.  L.  Veal,  "The  Dimension  of  Ice  Crystals  In 
Natural  Clouds,"  J.  Atmos.  Scl.  27,  919  (1970). 

32.  K.  0.  L.  F.  Jayaweera  and  B.  F.  Ryan,  "Terminal  Velocities  of  Ice 
Crystals,"  Quart.  J.  Roy.  Met.  Soc.  98,  193  (1972). 

33.  K.  0.  L.  F.  Jayaweera  and  B.  J.  Mason,  "The  Behavior  of  Freely  Fall¬ 
ing  Cylinders  and  Cones  In  a  Viscous  Fluid,"  J.  Fluid  Mech.  2£,  709 
(1965). 

34.  G.  M.  Bragg,  L.  van  Zuider,  and  C.  E.  Hermance,  "The  Free-Fall  of 
Cylinders  at  Intermediate  Reynolds  Numbers,"  Atmos.  Environ.  8,  755 
(1974). 

35.  K.  0.  L.  F.  Jayaweera  and  B.  J.  Mason,  "The  Falling  Motions  of  Loaded 
Cylinders  and  Discs  Simulating  Snow  Crystals,"  Quart.  J.  Roy.  Met. 

Soc.  92,  151  (1966). 

36.  J.  Zikmunda  and  G.  Vali,  "Fall  Patterns  and  Fall  Velocities  of  Rimed 
Ice  Crystals,"  J  Atmos.  Sci .  29 [,  133*  (1972). 

37.  R.  P.  Whitten,  "An  Investigation  of  Some  Aerodynamic  Factors  Affect¬ 
ing  Meteorological  Instrument  Readings  on  a  C130A  Research  Aircraft," 
Allied  Research  Associates,  AFRD  TN-60-454  (15  May  1960). 

38.  R,  D.  Ingebo,  "Drag  Coefficients  for  Droplets  and  Solid  Spheres  in 
Clouds  Accelerating  in  Alrstreams,"  NACA-TN  3762  (Sept.  1956). 

39.  T.  L.  Ogden  and  K.  0.  L.  F.  Jayaweera,  "Drag  Coefficients  of  Water 
Drops  Decelerating  in  Air,"  Quart.  J.  Toy.  Met.  Soc.  97,  571  (1971). 

40.  C.  N.  Davies,  "Definitive  Eguations  for  the  Fluid  Resistance  of 
Spheres,"  Proc.  Phys.  Soc.  (London)  J7,  259  (1945). 

41.  M.  Kajikawa,  "A  Model  Experimental  Study  on  the  Falling  Velocity  of 
Ice  Crystals,"  J.  Met.  Soc.  Japan  4£,  367  (1971). 

42.  K.  0.  L.  F.  Jayaweera  <ind  R.  E.  Cottis,  "Fall  Velocities  of  Plate-Like 
and  Columnar  Ice  Cryst.  Is,"  Quart.  J.  Roy.  Met.  Soc.  95,  703  (1969). 


-  101 


